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Abstract 

Orbital parameters of binary radio pulsars reveal the history of the pul- 
sars' formation and evolution including dynamic interactions with other 
objects. Advanced technology has enabled us to determine these or- 
bital parameters accurately in most of the cases. Determination of post- 
Keplerian parameters of double neutron star binaries (especially of the 
double pulsar system PSR J0737-3039) provide clean tests of the general 
theory of relativity and in the future may lead us to constrain the dense 
matter Equations of State. For binary pulsars with main-sequence or 
white dwarf companions, knowledge about the values of the orbital pa- 
rameters as well as of the spin periods and the masses of the pulsars and 
the companions might be useful to understand the evolutionary history 
of the systems. As accreting neutron star binaries lead to orbit circular- 
ization due to the tidal coupling during the accretion phase, their descen- 
dants i.e. binary millisecond pulsars are expected to be in circular orbits. 
On the other hand, dense stellar environments inside globular clusters 
cause different types of interactions of single stars with pulsar binaries 
like "fly-by", "exchange", and "merger". All these interactions impart 
high eccentricities to the pulsar binaries. So it is quite common to get 
eccentric millisecond pulsar binaries in globular clusters and we find that 
"fly-by" causes intermediate values of eccentricities while "exchange" or 
"merger" causes high values of eccentricities. We also show that "ioniza- 
tion" is not much effective in the present stage of globular clusters. Even 
in the absence of such kinds of stellar interactions, a millisecond pulsar 
can have an eccentric orbit as a result of Kozai resonance if the pulsar 
binary is a member of a hierarchical triple system. PSR J1903+0327 is 
the only one eccentric millisecond pulsar binary in the galactic disk where 
stellar interactions are negligible. The possibility of this system to be a 
member of a hierarchical triple system or past association of a globular 
cluster have been studied and found to be less likely. 
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1 Introduction 



Radio pulsars are rotating, magnetized neutron stars (NSs) emitting radio beams 
from magnetic poles. Their masses are in the range of 1 — 2 Mq, radii ~ 10 
km, densities ~ 10 14 g cm 3 and surface values of the dipolar magnetic field are 
around 10 s — 10 14 gauss. Presently there are 1864 radio pulsars with their spin 
periods covering a wide range 1.39 ms — 11.78 sQ. This number is increasing 
quite rapidly after the first serendipitous discovery in 1967 by Jocelyn Bell, a 
graduate student in Cambridge pQ. 

Pulsars arc believed to be born with spin period ~ 0.1 s which increases very 
slowly with time (at the rate of 10~ 20 — 10~ 14 ss~ 4 ) as their rotational energy 
is emitted in the form of electromagnetic energy. Millisecond spin periods are 
explained as a result of spin-up process, which occurs when a neutron star has 
a normal star as its binary companion and accretes matter from the companion 
(at its red giant phase) causing transfer of angular momentum. Depending upon 
the evolution of the companion, the resultant millisecond pulsar can be either 
isolated or a member of a binary system. In this article we concentrate on binary 
pulsars and discuss how the knowledge of their orbital parameters help us to 
understand the structure, formation, evolution and dynamic history of pulsars. 
As an example, binary millisecond pulsars are expected to be in circular orbits 
as tidal coupling during the accretion phase lead to orbit circularization. But 
if there is no history of accretion phase, the orbit of the binary pulsar can be 
eccentric because of the asymmetric kick imparted to the pulsar during its birth 
in a supernova. A millisecond pulsar binary can also have an eccentric orbit if 
it experiences interactions with other near-by stars (which is possible in dense 
stellar environments like globular clusters) or if it is a member of a hierarchical 
triple. 

It is noteworthy to mention here that in addition to radio wavelengths, pul- 
sars emit in x-rays and/or gamma rays too. But radio pulsars are the most 
useful astronomical objects firstly because the underlying neutron stars are 
more stable than in the case of X-ray pulsars which have episodically vary- 
ing accretion torques; secondly high sensitivities of radio telescopes and data 
analysis software make it possible to determine their spin periods very accu- 
rately. Moreover, for radio pulsars in binary systems (~ 8% of presently known 
radio pulsars are in binary systems) , it is possible to determine the mass of the 
pulsar and its companion, size and shape of the orbit etc by modelling Keplerian 
and post-Keplerian parameters. Such accurate determination of those parame- 
ters are difficult for X-ray pulsars at present. That is why we concentrate only 
on binary radio pulsars when we try to use the knowledge about their orbital 
parameters to understand the properties of pulsars. The first binary pulsar PSR 
B1913+16 was discovered in 1974 at Arecibo observatory by R. Hulse and J. H. 
Taylor 2 which is a 59 ms pulsar having orbital period 0.323 days and orbital 
eccentricities 0.671. Interestingly, its ~ 1.06 Mq mass companion star is most 
probably another neutron star i.e. the first discovered binary pulsar is also the 

1 http: / / www. atnf . csiro.au / research /pulsar /psreat/ 



2 



first discovered double neutron star (candidate) system. Afterwords nine other 
double neutron star (DNS) systems have been discovered. Among these ten 
DNSs, seven systems are confirmed to be DNSs and the rest are most likely 
to be DNSs [3]. These DNSs are useful tools to test the validity of the gen- 
eral theory of relativity (2] and even in principle to constrain the dense matter 
Equations of State [5]. In Section [5] we discuss how the study of DNSs can lead 
to a better understanding of the structure of neutron stars (i.e the Equation of 
State), the physics of the pulsar magnetosphere and the formation scenario for 
the special case of the double pulsar system. 

Most of the binary pulsars have either main sequence (MS) or white dwarf 
(WD) companions. Knowledge about the values of the orbital parameters as well 
as of the spin periods and the masses of these pulsars (and their companions) 
might be useful to understand the evolutionary history of these systems. The 
possible causes of eccentricities of millisecond pulsar binaries have been studied 
in this article. When these systems are in dense stellar environments as in 
globular clusters, the eccentricities can arise as a result of different types of 
interactions of single stars with pulsar binaries like "fly-by", "exchange" and 
"merger" (Section [3]). The eccentricity of the only eccentric millisecond pulsar 
PSR J1903 + 0327 in the galactic disk is difficult to explain as the hierarchical 
triple scenario faces difficulties to satisfy observed parameters of this system 
and the globular cluster origin also appears to be unlikely (section 2]) . All the 
findings of the present work have been summarized in section [5] 

2 Double Neutron Star Binaries 

After performing pulsar observations, one determines "time of arrivals" (TOAs) 
of pulses at the telescope. Then various clock corrections are applied to the 
TOAs, followed by some other corrections e.g. "Romer delay", "Shapiro delay" 
and "Einstein delay" to calculate the TOAs with respect to the solar system 
barycenter. Then the correction for the frequency dependent delay (dispersion) 
by the interstellar medium is added. Finally, the pulse phase <\> is fitted with 
time as a Taylor series expansion like <f> — tfio + v(t — to) + kv(t — to) 2 where 
4>o and to are reference phase and time, v and v are the spin frequency and 
its first time derivative. In case of a binary pulsar, one needs to model the 
orbital motion of the pulsar, and before doing that same type of delays for 
the binary system (e.g. another set of "Romer delay", "Shapiro delay" and 
"Einstein delay") should be taken into account. For a non-relativistic binary 
pulsar, five Keplerian parameters are required to describe the orbit, these are 
the orbital period (P or b), projected semi-major axis of the pulsar orbit (a p sin i, 
where i is the inclination of the orbit with respect to the sky plane), orbital 
eccentricity (e), longitude of periastron (w) and the epoch of the periastron 
passage (To). The semi major axes of the pulsar orbit and the companion orbit 
can be written as a p — am c /(m p + m c ) and a c = am p /(m p + m c ) respectively 
where a is the semi major axis of the relative orbit coming in the Kepler's law 
P^ rb /47r 2 = a 3 /G(m p + m c ). Here m p is the mass of the pulsar and m c is the 
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mass of the companion. Pulsar timing analysis means first getting an initial 
fit for the orbital parameters by inspecting the variation of the pulse period 
(which is equal to the spin period P s of the pulsar) and then improving the 
values of the parameters through a phase connected timing solution. But for 
relativistic binaries {i.e. when the companion of the pulsar is another neutron 
star or a white dwarf), five post-Keplerian (PK) parameters are needed, which 
are the relativistic advance of periastron (u> gr ), redshift parameter (7), Shapiro 
parameters (range "r" and shape "s" ) and the rate of change of orbital period 
(P or b) due to gravitational wave emission. Details of all the above mentioned 
parameters and how they are used in the timing models arc described very nicely 
in "Handbook of Pulsar Astronomy" by D. R. Lorimer and M. Kramer [BJ. 

Post-Keplerian descriptions of the orbits are essential for "double neutron 
star" (DNS) binaries which include PSR J0737-3039(A, B), PSR B1534+12, 
PSR J1756-2251, PSR J1811-173, PSR J1906+0746, PSR B1913+16, PSR 
B2127+11C (confirmed DNSs), PSR J1518+4904, PSR B1820-11, PSR J1829+2456 
(candidate DNSs - orbital parameters hint DNS nature but companion mass in- 
formations are insufficient to confirm). Only one of the DNSs, namely PSR 
J0737-3039 is a double pulsar system i.e. both the neutron stars emit radio 
pulses. This feature has made this system a unique astrophysical laboratory 
and we shall discuss about it in more details in subsections 12.11 12.21 and 12.31 
Moreover, there are three relativistic NS-WD systems - PSR J0751+1807, PSR 
J1757-5322, PSR J1141-6545 for which also post-Keplerian descriptions are un- 
avoidable. Let us now try to understand how DNSs can help us in better 
understanding of the structure of neutron stars, physics of the pulsar magneto- 
sphere and formation scenario (for the special case of the double pulsar system) 
of neutron stars. 

2.1 Constraining Dense Matter Equations of State 

After the first proposal of neutron star as the ultimate fate of stars (of mass 
~ 5 — 60 Mq) by Baade and Zwicky in 1934 [7] and the immediate prediction of 
mass and radius by Oppenheimer and Volkoff in 1939 [8], many people worked to 
understand the nature of matter at such high densities, resulting a large number 
of Equations of States (EsoS) including the EsoS for strange quark matter. 
Different EsoS result different mass-radius relations - see Fig 2 of [5] to realize 
how different the mass-radius curves can be. So the question arises, how to test 
the validity of those EsoS. Presently there are a number of efforts to constrain 
the EsoS through astronomical observations in X-rays. The usual approach 
is to determine the mass and the radius of the stars with the help of various 
observational features like gravitational redshifts (z) from spectral lines, cooling 
characteristics, kHz quasi-periodic oscillations (QPO) etc [51 [TU1 fTTT [T2l H"4"] . 
But these methods are not foolproof, e.g. the value of z used in Ozel's [TT] 
analysis of EXO 0748—676 could not be reproduced afterwords [15]. Moreover, 
to constrain EsoS from QPO observations, one needs to believe in a particular 
model of QPO which is again a subject of debate. An alternative method to 
constrain the dense matter EsoS might be the measurement of the moment of 
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inertia of a neutron star in a relativistic binary, especially for a DNS system and 
to match this observationally determined value with the theoretically predicted 
one. The Equation of state (EoS) which gives the best agreement between 
observed and the theoretical values of the moment of inertia of the NS can be 
considered as the best description of NS structure. 

The moment of inertia of a neutron star of known mass can be theoretically 
calculated for any particular EoS either using Hartle's approximations [16] or 
by performing exact calculations [171 118) for the rotating axisymmetric neutron 
star. It has been checked that the values of moment of inertia for a particular NS 
obtained by this two approaches do not differ more than 1% for a star like PSR 
J0737-3039A [T3Q1JEEE]. We found the moment of inertia of PSR J0737-3039A 
to be in the range 0.67 — 0.74 x 10 45 g cm 2 for a set of EsoS having different 
stiffness [T9] . 

Let us now discuss about observational determination of the moment of a 
neutron star by accurate measurement of the periastron advance (u! gr ) of the 
orbit of a relativistic binary pulsar (preferably a DNS) through timing analysis. 
The theoretical expression of ui gr including the effect of spin-orbit coupling is 
known to be as follow 1201: 



3/3 2 O or {, 
1 -e 2 



[1 + Ml - (ff.ij9.iA) + ff s2 j3 s2 /3o)] (1) 



where the first term inside the square bracket is the 1PN term, the second term 
is the 2PN term, and the third term represents the spin-orbit coupling (SO). 
The parameters used in the above expression are defined as : 
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where G is the gravitational constant, I a is the moment of inertia of the a th 
body (a = 1, 2), cu a = 2tt/ P S:a is its angular velocity of rotation (here a + 1 
means modulo 2, i.e. 2+1=1) and Q or b = 2ir/P or b. P s ^ a is the spin period of the 
a th star and P or b is the orbital period of the binary. x\ — nii/M, xi — m^jM 
and M — m\ + mi- k is the orbital angular momentum vector, s a is the spin 
vector, Kq is the line of sight vector perpendicular to the sky and i is the 
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inclination of the orbit with respect to the sky implying that the angle between 
Ko and k is also i. 

In the expression of Co gr , the moment of inertia of both members of the DNS 
system are present (through j3 Si i and f3 s , 2 in Eqn. [T]) which means that we 
have only one observable and two unknowns making Eqn. (fTJ) unsolvable. But 
fortunately, (3 sa is negligible for P SiQ > 100 ms. So if P s , 2 > 100 ms, equation 
for periastron advance becomes 



3$ Q 



orb 



[1 + Ml - [g A PM] (6) 



So DNS binaries having one fast neutron star (P s < 100 ms) and another slow 
neutron star (P s > 100 ms) are needed to extract the moment of inertia of 
the faster one by equating the observed uj gr with the prediction of Eqn. ([6j 
using measured values of e, P or b, h mi, m 2 an d Ps,i- It is clear that the angle 
between Si and k plays a vital to in the expression of g s \. If Si || k, the angle 
between K and Si becomes i and we get g s iu = a:i ^ 1 ^/^ ■ But if this is 
not the case i.e. the angle between si and k has a non-zero valne, g s i will be 
maximum if Si is parallel to the vector [(3 sin 2 i — 1) k . si + cos i Ko . Si] 

1 /2 

giving | g shmax \= [3 + -X^} 



for other orientation of s 1; g s i 
So the knowledge of the angle 



3:1(43:1+33:2) 

[ Z{l-e 2 ) 1 /' 2 _ 

will be between + | g s i ima x | to — | g s i, ma x 
between k and s 1 is needed for determination of the moment of inertia of the 
faster member of a DNS system. 

At present, it seems that the use of Eqn. ^ to determine the moment of 
inertia of a neutron star is possible only for the double pulsar system PSR J0737- 
3039 [HI [19] , as for this system, timing analysis of the pulses from both the stars 
has already enabled us to determine e, P or b, mi, m%, P Sj i and P s ^ [22j|23]. The 
two pulsars which are named as PSR J0737-3039A and PSR J0737-3039B have 
masses and spin periods as tua — 1.3381 M©, tub — 1.2489 M , P s _a = 22.7 
ms and P S) b = 2.77 s with orbital parameters as e = 0.0877, P or b — 0.102252 
days, P or t = —1.252 x 10~ 12 , i — 88.69°. So it is clear that we can easily 
neglect the spin-orbit term due to pulsar B. Although the angles between the 
vectors k, Si and Ko, Si are not known which are needed to be determined 
first to measure the moment of inertia of the pulsar A accurately. Just for the 
sake of completeness, we wish to mention that for PSR J0737-3039, the value 
of I g s A,max | is 1.2176 whereas the value of g s A,\\ is 0.60876. The accuracies of 
other relevant parameters are also needed to be improved. 

In general, it is clear that the {gsiPsiPa) term in Eqn. ^ will have higher 
value for (i) smaller value of P s ,i, (ii) smaller value of P or b and (iii) higher value 
of e. Unfortunately, the eccentricity of PSR J0737-3039 is small in comparison to 
other DNSs. Future discovery of any double pulsar system satisfying these three 
conditions better than PSR J0737-3039 will be useful for the purpose of moment 
of inertia measurement (in addition to the necessary condition P Sj 2 > 100 ms). 



6 



2.2 Interplay between Two Members of the Double Pulsar 



Among many interesting features observed from the double pulsar system PSR 
J0737-3039, one is the variation of the flux density of the slower pulsar B at 
its different orbital phases due to the distortion of B's magnetosphere by A's 
wind [231 [53]. On the other hand, the flux density of pulsar A is constant 
for most of its orbit except 30 s near its superior conjuction when the radio 
beam from the pulsar A is eclipsed by the magnetosphere of the pulsar B [25] . 
The eclipse is asymmetric, with the flux density decreasing more slowly during 
ingress than it increases during egress. The duration of the eclipse is very mildly 
dependent on the observing frequencies, with the eclipse lasting slightly longer 
at lower frequencies [25] ■ Theoretical models explain these eclipse properties 
in the context of synchrotron absorption of A's radio emission by the plasma 
in the closed field line region of B's magnetosphere [23 [23 [2E]- By modeling 
the eclipse, the precession of B's spin axis about the orbital momentum vector 
has been found to match with the prediction of the general theory of relativity 
within the 13% measurement uncertainty |29j . It is noteworthy to mention here 
that till date, the eclipse has been studied in the observing frequency range 
427-1400 MHz. The study of PSR J0737-3039 (both A and B) at a much lower 
frequency of 326 MHz is being performed using Ooty Radio telescope (RAC- 
TIFR, Ooty, India). The details of eclipse characteristic at this low frequency 
will be useful to test the validity of the eclipse models and better understanding 
of the pulsar magnetosphere. 

2.3 Formation of the Double Pulsar 

Using observed parameters of the double pulsar system, different formation 
models were proposed [30l [31] [32] . Initially, it was proposed that this system 
has originated from a close helium star plus a neutron star binary in which at 
the onset of the evolution the helium star had a mass in the range 4.0 — 6.5 
Mq and an orbital period in the range 0.1 — 0.2 day [30]. In this model, a kick 
velocity in the range 70 — 230 km s _1 must have been imparted to the second 
neutron star at its birth. Afterwords an electron capture supernova scenario 
was proposed [31]. At the same time, Piran & Shaviv [32] estimated that it is 
most likely that the progenitor of the pulsar B had a mass of 1.55 ± 0.2 M Q 
with a low kick velocity 30 km s _1 . On the other hand, a study based on the 
population synthesis method [33| favored the standard formation scenario with 
reasonably high kick velocities (70-180 km s _1 ). Later, it has been suggested 
[19] that the pulsar B is a strange star rather than a neutron star as the baryon 
loss is smaller for the first case supporting the models with small kick. 

3 Binary Pulsars in Globular Clusters 

Globular clusters (GCs) are dense spherical collection of stars orbiting around 
the galactic center containing a significant number of binary stars. The first 
radio pulsar discovered in a globular cluster is PSR B1821-24A in M28 34 in 
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1987, 20 years after the discovery of first pulsar by Jocelyn Bell. This pulsar is 
at present an isolated pulsar but its very low spin period (3 millisecond) hints to 
a binary origin and spin-up due to accretion . How old and evolved stellar sys- 
tems like GCs can contain apparently young, active objects like pulsars have led 
to suggestions, based on lifetime constraints from orbital decay by gravitational 
radiation [35]; that some of these systems have been formed recently, or are 
forming even today. Interaction between binary stars and single stars in GCs is 
believed to be an important dynamical process as the resulting binary systems 
may provide a substantial source of energy for the host GC, as the binding en- 
ergies of a few, very close binaries (like neutron star binaries) can approach that 
of a moderately massive GC ISZ] • Observable parameters of the pulsars and 
their binary companions in GCs, such as spin, orbital period and eccentricity, 
projected radial position in the cluster, companion mass and their distributions 
provide a tracer of the past history of dynamical interactions of the binary NSs 
in individual GCs. These parameters can even provide a valuable test-bed to 
examine the theoretical scenarios of formation and evolution of recycled pulsars. 

3.1 Archival Data of Binary Pulsars in Globular Clusters 

Till today, a total of 140 pulsars in 26 GCiS have been discovered. Among them 
74 are known as binaries (in 23 GCs), 59 are known as isolated and 7 have no 
published timing/orbital solutions. Orbital parameters of one binary, namely 
PSR J2140-2310B in M30 are not well determined, only lower limits of P orb 
and e are known as P or h > 0.8 days and e > 0.52. So we exclude it from our 
analysis and proceed with the rest 73 binaries. 

In Fig [TJ we plot e Vs P or b for these 73 binary pulsars in GCs with known 
orbital solutions. A logarithmic scale in both eccentricity and orbital periods 
are chosen, as the enormous range of both variables and the regions occupied 
by observed pulsars are less obvious in linear scales. Downward arrows have 
been drawn for the systems for which only upper limits of eccentricities (see 
table [2]) are known. Non-millisecond pulsars are enclosed by s an d the DNS 
system PSR B2127+11C in M15 is enclosed by a green □. In the present work, 
we define millisecond pulsars as the pulsars having P s < 30 ms, we don't go for 
other definitions of millisecond pulsars based on both P s and P s (see [38] for one 
such definition) as for GC binaries, the values of P s are effected significantly by 
the cluster potential. 

Observed binaries in GCs can be categorized in three groups : (I) large 
eccentricity pulsars (e ^ 0.01), 21 pulsars (22 if we include M30B), (II) moderate 
eccentricity pulsars (0.01 > e 3 x 10~ 6 ), 20 pulsars and (III) small eccentricity 
pulsars (e ~ 0), 32 pulsars. In the database, orbital eccentricities of group (III) 
pulsars have been listed as zero, and in this logarithmic plot we assign them an 
arbitrarily small value ofe = 3 x 10~ 7 . 

In addition we plot on the right of Fig [U dashed lines corresponding to 
the freeze-out eccentricity - orbital period relation predicted by Phinney |39j 

information on these pulsars are found from P. Freire's webpage updated August 2008, 
http: / /www. naic.edu/~pfreire/GCpsr. html 
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Fi gure 1: Plot of eccentricity vs orbital periods for binary radio pulsars in globular clusters. 
The DNS system PSR B2127+11C ib M15 is enclosed by a green □. Non-millisecond pulsars 
are enclosed by O s - The dashed lines in the right side represents the eccentricity-orbital 
period relation as predicted by Phinney |39| - the upper line is the original line drawn by 
Phinney 39_ and the lower line is using population II models. The black dotted lines parallel 
to the x-axis have been drawn to distinguish three groups of GC pulsars. 

on the basis of the fluctuation dissipation theorem for convective eddies in the 
erstwhile red giant envelopes surrounding the white-dwarf cores which end up 
as companions of the neutron stars. This relation is obtained from the following 
expression: 



where Vtorb = 2-Tr / P or b is the orbital angular rotation frequency, a is the semi- 
major axis of the binary, e is the orbital eccentricity, L is the luminosity of 
the companion, R c is the radius of the companion, m env is the envelop mass 
of the companion, m p is the pulsar mass, m c is the companion mass and ji = 
(m p m c )/(mp + m c ) is the reduced mass of the system. Phinney used R c — 
0.125 a when the mass transfer ceases and e freezes. Then the above equation 
reduces to < e 2 >V 2 = C x (P or6 /100) where C x is a function of the effective 
temperature (T e j^), [i and m env . P or h is in days. Phinney calculated C\ — 1.5x 
10~ 4 (sec Eqn. 7.35 of [35]) using average values of m env and T e // f° r population 
I stars (z = 0.04) [40] which gives the upper dashed green line in the right side 
of Fig[T] However for globular clusters, population II stars are more appropriate. 
So we recalculated C\ for population II red giant stars [41] and obtained C\ = 




(7) 
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0.86 x 10~ 4 . This line is plotted as the lower dashed green line in Fig. [TJ 

Earlier, it had been found that while the low mass X-ray binaries (LMXBs) 
are predominantly found in very dense GCs, the radio pulsars tend to be found 
more evenly distributed among low and high density clusters 42 . We also don't 
find any obvious correlation of the number of binary radio pulsars with any GC 
parameters like d, i>io, n± or Vio/n^ (where 714 is the number density (n) of 
single stars in units of 10 4 pc -3 and v\q is the velocity dispersion (v) in units of 
10 km s _1 in GCs and d is the distance of the GCs from sun in kpc see Table 
[1]). But if we take only GCs with more than 3 observed binary pulsars, we find 
the number of binary pulsars to increase with n^, but there are only five such 
GCs (Table [J). 

3.1.1 A comparison of binary pulsars in globular clusters and binary 
pulsars in the galactic disk 

There are total 83 binary pulsars in the galactic disk (and 1652 isolated pulsars) 
among which orbital period and/or eccentricity values are unknown for four 
binaries. So we exclude them from our study and we consider the rest 78 binary 
pulsars. In the database, eccentricities of three binaries have assigned the value 
"0" for which we put an arbitrary small value 3 x 10 -7 (as in the case of GC 
binary pulsars). 

In Fig[21 we plot e Vs P or j, for all binary pulsars with known orbital solutions. 
Blue Ts are for GC binaries and red As are for disk binaries. Three disk binaries 
with P OT b > 1000 days have not been plotted here as we have selected the range 
of P or b up to 1000 days in this plot. Here we take upper values of eccentricities 
as the actual values. DNSs are enclosed by green Ds. Non-millisecond pulsars 
are enclosed by s - It is believed that millisecond pulsars in binary systems 
formed from mass and angular momentum transfer due to Roche lobe overflow 
and the resultant tidal effects should appear in e ~ 10~ 6 — 10~ 3 [3S]. Since many 
highly eccentric binary millisecond pulsars are found in GCs, this indicates that 
stellar interactions are important for inducing higher eccentricities which will 
be discussed later in this section. On the other hand, PSR J1903 + 0327 is the 
only one eccentric millisecond pulsar in the disk and there are lots of works in 
the recent past to explain the eccentricity of this system [33 H31 ES] ■ We shall 
discuss about this system in somewhat details in section [4] 

Observational selection effects may be influencing the distribution seen in 
this figure for the two sets [42], the most important selection effect being op- 
erative towards the left of the diagram, namely, it is more difficult to detect 
millisecond pulsars in very short orbital period and highly eccentric binaries; 
distance to the pulsars also is an important selection effect, since at large dis- 
tance only the brightest pulsars can be observed. Nevertheless, it is clear that 
there is a large abundance of pulsars with long orbital periods and intermedi- 
ate eccentricity near the eddy "fluctuation dissipation" lines among the galactic 
disk pulsars which are absent in the GC pulsar population, where the important 
selection effects are unlikely to play a major role. The observed preponderance 
of DNSs in the disk with the very short orbital periods and large eccentricities 
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Figure 2: Plot of eccentricity vs orbital periods for binary radio pulsars in and out globular 
clusters. Blue Ts are for globular cluster binaries and red As are for disk binaries. Double 
neutron star binaries are enclosed by green CDs. Non-millisecond pulsars are enclosed by O s - 
Three disk binaries with P or t > 1000 days have been excluded. The dashed lines in the right 
side represents the eccentricity-orbital period relation as predicted by Phinney [39| - the upper 
line is the original line drawn by Phinney 1391 and the lower line is using population II models. 
The black dotted lines parallel to the x-axis have been drawn to distinguish three groups of 
GC pulsars. 

may be related to selection effects, but there is at least one DNS in GC - PSR 
J2127+1105C in M15 in the same phase-space. The disk DNSs are by and large 
at small distances while M15 and other host GCs are usually far away. 

In the top panel of Fig. [3] we compare the distribution of m c s for the GC 
and the disk binaries and found that there are more low mass companions for 
GC binaries. In the bottom panel of Fig. [3] we compare the distribution of P s s 
for the GC and the disk binaries and found that in GCs, most are millisecond 
pulsars while for disk pulsars, P s are more uniformly distributed (not completely, 
more pulsars with lower P s ) over a wide range. 
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Figure 3: The distribution of companion masses for the binary radio pulsars (top) and the 
distribution of spin periods for the binary radio pulsars (bottom). Left panels are for the GC 
binaries and the right panels are for the disk binaries. 
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In the left panel of Fig. @J we plot P or b, e and m c for all binaries while in 
the right panel we plot P or b, e and P s . In both cases, blue points represent 
GC binaries and red points represent disk binaries. Three disk binaries with 
m c > 2.0 Af (see tabled have been omitted from the plot in the left panel, even 
then we can see that companion masses of disk binaries are usually higher than 
those of GC binaries. Nearly circular binary pulsars have low mass companions. 




Figure 4: Left panel shows the variation of companion masses with P or i, and e, blue points 
are for GC binaries and red points are for disk binaries omitting three disk binaries with 
m c > 2.0 MQ(see table [3} for the ease of comparison. Right panel shows the variation of P a 
with P or i and e, blue points are for GC binaries and red points are for disk binaries. Note 
the presence of eccentric millisecond binaries in GCs. 

It was proposed that orbital parameters of binary pulsars with white-dwarf 
companions satisfy the following relation |46j : 

Porb (1 - ef' 2 = 0.374 

In the upper panel of Fig. El we plot P or b (1 — e) 3 / 2 against m c for binary 
pulsars in GCs and also in the galactic disk. DNS binaries and the binaries 
with m c > 1.44 M Q (Chandrasekhar limit) have been excluded. GC binaries 
are shown with blue Ts and disk binaries are shown with red As. The points 
correspond to median values of m c taking i = 60° while the upper and the 
lower limits of m c s correspond to i = 26° and i — 90° respectively. Theoretical 
prediction following Eqn. JHJ is the middle solid green line. The upper and lower 
lines are obtained by multiplying the right hand side of Eqn. (JSJ by 2.5 (see [4"B]). 
Note that the model used to derive Eqn. (|SJ) is not valid when m w d < 0.13 M Q , 
though we have plotted the binaries with m c < 0.13 M (located left of the 
vertical black dotted line) just to see how P or b (1 — e) 3 / 2 is related with m c for 
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these systems. For recycled (millisecond) pulsars, where significant amount of 
mass had been transferred to the pulsars in past, one should expect lower values 
of m c and e ~ 10 -6 — 10~ 3 . for these cases the relation becomes 



0.374 



4.5 



4950 ml 
1 + 4 mi 



0.5 



3/2 



m . 



-1/2 



(9) 



The deviation of GC binary pulsars from the theoretical prediction of Eqn. 
ijHJl is well understood as their eccentricities have been effected by stellar inter- 
actions. So for the millisecond pulsar binaries in GCs (which naturally exclude 
the DNS in M15), we plot Eqn. ((9]) in the lower panel of Fig. (J5j) . Even then 
we found significant discrepancy between the theoretical prediction and the ob- 
served values. Moreover, for binary pulsars in the disk (upper panel of Fig. [5j 
where eccentricities and orbital periods should bear the clear signature of stellar 
evolutions, we also found disagreement between the model and the observation. 

We summarize the relevant properties of the host GCs and in Table (fT|). 
The properties of binary pulsars in globular clusters and in the galactic disk are 
given in tables (|2|) and ([3]) respectively. Interested readers are requested to see 
original sources for better precision of the given parameters and for the values 
of other parameters. 

Through a combination of two-body tidal captures and exchange interactions 
where an isolated neutron star replaces a normal star in a binary, numerous 
neutron star binaries are formed in GCs |47j . This explains why the ratio of 
binary pulsars to isolated pulsars is so high in GCs in comparison to that in the 
galactic disk. Tidal capture leads to tight binaries whereas exchange usually 
leads to wider binaries (see Eqns. [T4|) . If the mass transfer in the tidally 
captured neutron star binary happen in a stable manner, it can spin up the 
neutron star to a millisecond spin period. Thus the observed preponderance of 
millisecond pulsars in GC binaries can be attributed to the tidal capture process 
(details will be discussed in the next subsection). We explain the eccentricities 
of millisecond pulsar binaries in GCs as a result of different types of interactions 
of single stars with pulsar binaries using present properties of the GCs. Some of 
the results can be found in our earlier papers [2H1 HI] • For the sake of simplicity, 
we don't go for a full population synthesis study of the evolution of the globular- 
clusters including all kinds of stellar interactions e.g. binary-binary, binary- 
single star and single star - single star interactions. We consider only binary- 
single star interactions faced by binary pulsars preceded by a discussion on the 
formation of binary pulsars through tidal capture. 
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Figure 5: Upper Panel : Plot of P orb (1 - e) 3 / 2 against m c for binary pulsars in GCs (blue 
Ts) and in the disk (red As). Lower Panel : Ifflot of P or t against m c for binary millisecond 
pulsars in GCs. In both of the panels, DNSs and binaries with m c > 1.44 Mq have been 
excluded. The points correspond to median values of m c taking i = 60° while the upper 
and the lower limits correspond to i = 26° and i = 90° respectively. Theoretical prediction 
following Eqn. J5J or Eqn. (O are the middle lines. The upper and lower lines are obtained 
by multiplying the right hand side of the Eqns. (|8j [9j by 2.5 (see [46] ) . Note that the model 
used to derive the equations is not valid when m c < 0.13 Mq, though we have plotted the 
binaries with m c < 0.13 Mq (located left of the vertical black dotted line) in both the panels. 



3.2 Orbital Parameter Changes Due to Stellar Interac- 
tions and Gravitational Radiation 

Globular clusters contain a large number of LMXBs compared to the galactic 
disk; this had led to the suggestion 50 that a binary is formed by the tidal 
capture of a non-compact star by a neutron star in the dense stellar environment 
of the GC cores. If mass and angular momentum transfer from the companion 
to the neutron star took place in a stable manner, this could lead to recycled 
pulsars in binaries or as single millisecond pulsars [511 [52] [53] . However large 
energies may be deposited in tides after the tidal capture of a neutron star 
by a low mass MS star in a GC. The resultant structural readjustments of 
the star in response to the dissipation of the modes could be very significant 
in stars with either convective or radiative damping zones [541 155) and the 
companion star can undergo a size "inflation" due to its high tidal luminosity 
which may exceed that induced by nuclear reactions in the core, by a large factor. 
Efficiency of viscous dissipation and orbit evolution is crucial to the subsequent 
evolution of the system as viscosity regulates the growth of oscillations and the 
degree to which the extended star is bloated and shed. The evolution of the 
system could lead to either a merger (leading to a Thorne Zytkow object) or 
a neutron star surrounded by a massive disk comprising of the stellar debris 
of its erstwhile companion or, in other cases, an LMXB or even a detached 
binary following envelope ejection. Tidal oscillation energy can return to the 
orbit (i.e. the tide orbit coupling is a dynamical effect) thus affecting the orbit 
evolution or the extent of dissipative heating in the less dense companion. The 
evolutionary transition from initial tidal dynamics to any likely final, quiescent 
binary system is thus regulated by viscosity. Tide orbit coupling can even 
lead to chaotic evolution of the orbit in some cases [56] but in the presence of 
dissipation with non-liner damping, the chaotic phase may last only for about 
200 yrs after which the binary may undergo a periodic circularization and after 
about 5 Myrs finish circularizing [57]. Of the stars that do not directly lead 
to mergers, a roughly equal fraction of the encounters lead to binaries that 
either become unbound as a result of de-excitation or heating from other stars 
in the vicinity or they are scattered into orbits with large pericenters (compared 
to the size of the non-compact star) due to angular momentum transfer from 
other stars [58] . Thus "intermediate pericenter" encounters can lead to widened 
orbit (even 100 days orbital period in an eccentric orbit) in tidal encounters 
involving main sequence stars and neutron stars, which in the standard scenario 
are attributed to encounters between giants and neutron stars. The complexities 
of the dynamics of tidal capture and subsequent evolution are manifold and 
attempts to numerically simulate collisions of neutron stars with red giants have 
been made with 3D Smoothed Particle Hydrodynamics (SPH) codes [59l [60] . 
We note in this context that Podsiadlowski et al. [UJ have stated "it is not only 
premature to rule out tidal capture as a formation scenario for LMXBs, but 
that LMXBs in globular clusters with well determined orbital periods actually 
provide observational evidence in its favor". The different formation channels 
of pulsars in binaries in GCs and the birthrate problems of millisecond pulsar 
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binaries vs LMXBs were considered in [5TJ • 

Until the early eighties there was no evidence of a substantial population of 
primordial binaries in any globular clusters, and it was even thought that GCs 
are significantly deficient in binaries compared to a younger galactic population 
[55] . Theoretical modeling of GCs was often started off as if all stars were 
singles. During the 1980s several observational techniques began to yield a rich 
population of binary stars in GCs [64l [65l [66] . When the local binary fraction 
is substantial, the single star - binary interaction can exceed the encounter 
rate between single stars by a large factor [57]. The existence of a significant 
population of primordial binary population in GCs indicated that three body 
processes have to be accounted for in any dynamical study of binaries involving 
compact stars. An encounter between an isolated star and a binary may lead to 
a change of state of the latter, e.g.: i) the original binary may undergo a change 
of eccentricity and orbital period but otherwise remain intact - a "preservation" 
or a "fly-by" process; ii) a member of the binary may be exchanged with the 
incoming field star, forming a new binary - an "exchange" process; iii) two of 
the stars may collide and merge into a single object, and may or may not remain 
bound to the third star - a "merger" process; iv) three of the stars may collide 
and merge into a single object - a "triple merger" process; v) all three stars may 
become unbound - an "ionization" process. The ionization is always a "prompt" 
process, whereas the others can be either prompt or "resonant" process. In 
resonant processes, the stars form a temporarily bound triple system, with a 
negative total energy but which decays into an escaping star and a binary, 
typically after 10-100 crossing times. 

The value of the binding energy of the binary and the velocity of the in- 
coming star determine the type of interaction the binary will encounter. The 
critical velocity v c of the incoming star, for which the total energy (kinetic plus 
potential) of the three body system is zero, can be defined as 



, mim 2 (rai + m 2 + mz) 1 



vt = G v — ; — (io) 

m 3 (mi + m2) a in 

where m\ and mi are the masses of the binary members, m-$ is the mass of the 
incoming star, ai n is the semi-major axis of the binary and G is the gravitational 
constant. 

In case of a binary-single star interaction, semi major axis of the final binary 
(a fi n ) is related to the semi major axis of the initial binary as follows [57] 



1 m a m b 

a fin — i a ain V iJ J 

1 — A TO1TO2 

where m± and mi are masses of the members of the initial binary, m a and nib 
are masses of the members of the final binary. A is the fractional change of 
binary binding energy i.e. A = (Ei n — Efi n ) / Ei n . The binding energies of the 
initial and final binaries are 

_ ^ mmii _ mamb . , 

&in — > &fin — ^^TJ \ LZ ) 
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For fly- by, m a — mi, mj, = ni2 giving 

1 

a-fin — r _ ^ a m (13) 

For exchange, m a — mi, — 7773 (star 2 is being replaced by star 3) giving 
- 1 mim 3 _ 1 m 3 

1 — A miW2 1 — A 7712 

For merger, m a = mi, to?, = 777,2 + "13 (star 2 is being merged with star 3) giving 

1 7771(7712 + m 3 ) _L 7772 + 7773 , , 

Ifjn — - 7 CLin — ~ 7 din 

1 — A 77717772 1 — A 7772 

The actual value of A is not known. Putting A — simplifies equations (fTTj) . 
<jT3j) , (fT4|) and (fT5|) giving a/ in = s in for fly- by interactions. Similar expressions 
can be derived for other cases e.g. star 1 is being replaced by star 3 or star 1 
being merged with star 3. But in the present study, we always assume star 1 
to be the neutron star, so these processes are irrelevant. Merger interactions 
which do not result in binary systems are also irrelevant for this work. 

Interactions between binaries and single stars can enhance the eccentricity of 
the binary orbit and this will be discussed in detail in the subsequent sections. 
On the other hand, binary pulsars emit gravitational waves which reduce both 
the sizes and eccentricities of the orbit leading to mergers of the binary members 
on a timescale of t gw which can be calculated as follows [5S] : 



where 



and 



tgw = ~. — (16) 

&gw 

304 G 3 m p m c (rn p + m c 
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-9(e) (17) 



5(e)^(l- e 2 )" 5/2 e(l + ^e 2 j (18) 
which is almost same as the expression of t gw calculated from a gw . 

t gw = (19) 



and 



a 

64 G a 777p777 c (777p + 777 c ) 



gw 

" ' '/(e) (20) 



In this work, we use eq. (Tl6|) to calculate t 
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3.2.1 Fly-by interactions 



Figs. [TJ[2]and Table [2] show that most of the binary pulsars in GCs are millisec- 
ond pulsars (69 pulsars out of 73 have P s < 30 ms, one has P s = 33.16 ms and 
the other three having P s values as 80.34 ms, 111.61 ms and 1 s). Theoretically 
one expects that spun-up, millisecond pulsars in binary systems formed from 
mass and angular momentum transfer due to Roche lobe overflow and the resul- 
tant tidal effects should appear in low eccentricity orbits e ~ 10~ 6 — 10~ 3 [39] . 
Since many highly eccentric binary millisecond pulsars are found in globular 
clusters (see the right panel of Fig [4]), this indicates that stellar interactions are 
important in GCs for inducing higher eccentricities. Rasio and Heggie [551 [70] 
studied the change of orbital eccentricity (Se = e/j„ — Ri n where ei n is the initial 
eccentricity and e /,„ is the final eccentricity which is observable at present) of an 
initially circular binary (ej n = 0) following a distant encounter with a third star 
in a parabolic orbit. They used secular perturbation theory, i.e. averaging over 
the orbital motion of the binary for sufficiently large values of the pericenter 
distance r p where the encounter is quasi-adiabatic and used non secular per- 
turbation theory for smaller values of r p where the encounter is non-adiabatic. 
In the first case Se varies as a power law with r p /ai n and in the second case 
Se varies exponentially with r p /ai n . The power law dominates for eji n < 0.01 
wheras the exponential law dominates for e/, n > 0.01. The timescales for these 
processes are [69) (where we write e in place of e fi n for the sake of simplicity) : 



where P or b is in days and tfi y is in years. The values of via/n^ arc different for 
different GCs (table [TJ and we grouped them into six different groups accord- 
ing to their values of vio / 71,4 and calculated t fi y for mean values of vio /n^ for 
each group. The unit of is 10~ 3 km s _1 pc 3 which we will not mention 

explicitly anymore. 



t f iy = 4 x lO 11 ^ Vo^V/ 6 (e < 0.01) 
t flv = 2 x W^n^VwP-^ 3 Hn(e/4)]- 2/3 (e > 0.01) 



(22a) 
(22b) 
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Fi gure 6i Cluster binary pulsars in the e — P rb plane with contours of tfiy — 10*^ yrs (solid 
lines) and tji y = 10 s yrs (dashed lines) for each group. Contours of t gw = 10 10 yrs (solid lines) 
and t gw = 10 s yrs (dashed lines) for binaries with m p = 1.4 Mq and m c = 0.35 Mq or 
0.16 Mq are also shown. Pulsars with projected positions inside the cluster core are marked 
with +, those outside the cluster core are marked with X and the pulsars with unknown 
position offsets are marked with o. Solid red line corresponding to tfiy = 10 10 years for 
flo/ n 4 = 0.0024 is outside the range plotted. Pulsars in a particular group are marked with 
the same color as of the tfi y contours for thar group. The upper limits of eccentricities are 
taken as the actual values of eccentricities. 



Cluster binary pulsars in the e — P or b plane with contours of tfi y — 10 10 
yrs (solid lines) and tfi y = 10 s yrs (dashed lines) for each group are shown 
in Fig [5] Contours of t gw — 10 10 yrs (solid lines) and t gw — 10 8 yrs (dashed 
lines) for binaries with m p = 1.4 Mq and m c = 0.35 M© or 0.16 Mq are 
also shown. Pulsars with projected positions inside the cluster cores are marked 
with +, those outside the cluster cores are marked with x and the pulsars 
with unknown position offsets are marked with o. Solid red line corresponding 
to tfiy — 10 10 years for = 0.0024 is outside the range of axes used in 

our plots (away from the upper left corner). Pulsars in a particular group are 
marked with the same color as of the tfi y contours for that group. The upper 
limits of eccentricities are taken as the actual values of eccentricities. If a pulsar 
is located on the upper left half of the corresponding tfi y — 10 10 years line i.e. 
in the region where tfi y > 10 10 years, then its eccentricity cannot be due to 
fly-by interactions. There are two such binaries, one is PSRB1718-1^|( in NGC 
6342, in group 4) and the other is PSR B 1639+36B (in M13, in group 6). But 
the first one is a normal pulsar with P s = 1 s and as it is not spun up and 
may not have been subject to a great deal of tidal forces, it can have escaped 
circularization. On the other hand, the second pulsar is a millisecond pulsar 
with P s — 3.528 ms, so it should have been in a circular orbit if it has not had 
the time to go through any significant kind of stellar interactions. Therefore it 
is of interest to investigate if it is a result of an exchange or merger interaction. 
Moreover, some of the eccentric pulsars which seem to be explainable by fly-by 
interactions lie outside of globular cluster cores (where the stellar density is 
comparatively low and hence fly- by interactions are less efficient). For a few 
others positional offsets are not known. Even if a pulsar appears to be inside 
the cluster core in the projected image, it can still be outside the cluster core in 
the three dimensional space. GC models show that values of both vio and 
fall outwards from the cluster center and as the fall of 714 is much more rapid, 
the value of K10/714 is higher outside the cluster core. As we have already seen 
(different panels in Fig. [5]) that an increase in the value of ^10/^4 shifts the 
contour of tfi y = 10 10 years rightwards, this will make the pulsar fall in the 
region where tfi y > 10 10 years excluding the possibility of fly-by interaction 
as the cause of its eccentricity. In those cases also, we need to think about 
exchange and/or merger interactions. 

We have seen that there are many zero eccentricity pulsars which lie in the 
region where tfi y < 10 10 years. So it is of interest to understand why they 
still appear in circular orbits despite the possibility of fly-by interactions. We 
will discuss about this point in subsection 13.41 

3.2.2 Exchange and merger interactions 

Since exchange and merger interactions have not been amenable to analytic 
treatment one needs to use numerical techniques. We performed numerical 

3 Note however that its cluster association is sometimes doubted _42_, as it has a large offset 
from the GC center. 
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scattering experiments using the STARLAB task "sigma3" which gives the 
scaled cross section X for different types of interactions, e.g. fly-by, exchange, 
merger and ionization (both resonant and non- resonant). The inputs given are : 
velocity (v) of the incoming star (with mass and radius 7713, R3) and the masses, 
the radii, the semi major axis of the initial binary (mi, R\\ m^, R2] o-in) as well 
as the trial density which determines the number of scattering experiments to 
be done. X is defined as : 

AT = ^V(-V (23) 
7T af n \v c J 

where a is the cross section and v c is the critical velocity as defined earlier 
(eq. [TUJ) - So from the reported values of X, one can calculate a and then the 
interaction time scale ti„ t as U n t = 1/nav where n is the number density of 
single stars. We took different sets of stellar parameters (mi,Ri; m^Ri) to 
understand the effect of each of them and for each set of stellar parameters, 
we varied aj„ over a wide range in such a way that the initial orbital periods 
(Porb, in, calculated using Kepler's law) are in the range of 0.01 — 1000 days. 

For each set of parameters (mi , R\ ; mi,Ri; "in ) > we take the maximum trial 
density (n tr iai) as 5000 which is uniformly distributed in the impact parameter 
(p) over the range < p < po where po simply corresponds to a periastron 
separation of 2a. The impact parameter range is then systematically expanded 
to cover successive annuli of outer radii p, = 2 l / 2 p with n tr iai trials each until 
no interesting interaction takes place in the outermost zone [71] . Thus the total 
number of trials become significantly large. As an example, for the parame- 
ters mi = 1.4 M Q ,R l = 10 km; m 2 = 0.16 M e ,R 2 = 0.16 R G ; m 3 = 
0.33 M Q ,i?3 = Rq; a,i n = 0.2 AU; v = 11.76 km s _1 we had a sample total 
of 15570 scatterings. Out of these 10960 were fly-bys, 3722 exchanges, 982 two 
mergers and 6 three-mergers. 

We performed simulations over widely different values of v e.g. 11.76 km s _1 , 
13.15 km s _1 , 7.79 km s _1 and 3.29 km s^ 1 chosen in such a way that almost the 
entire range of U10/714 (0.0024, 0.062, 0.131, 0.934 respectively) has been covered. 
The variation in v is also significant. But we mainly concentrate on the runs for 
v = 11.76 km s _1 as in Ter 5, the host of the largest number of known binary 
radio pulsars. 

In Fig. {2} we plot the variation of X for exchange (+), merger (A) and ion- 
ization (o, whenever significant) with P or b, m and check the dependencies on dif- 
ferent parameters - (i) In the first figure, we set mi = 1.4 M , — 0.16 M , 
m 3 = 0.33 M Q , and take two values of v e.g. 13.16 km s _1 (red) and 3.27 
km s _1 (green) which are the highest and the lowest values of v among our choice 
(see Table [J), (ii) In the second figure, we set mi = 1.4 Mq, to 3 = 0.33 M Q , 
v = 11.76 km s _1 , and take different values of rri2 e.g. 0.40 Mq (red), 0.16 M Q 
(green) and 0.024 Mq (blue). Ionization starts for 7712 = 0.024 Mq at P or b, in ^ 
200 days, (iii) In the third figure, we set mi = 1.4 Mq, m 2 = 0.16 Mq, 
v = 11.76 km s _1 , and take different values of m 3 e.g. 0.33 Mq (red), 0.80 Mq 

4 www. ids. ias.edu/~star lab/ 
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(green), (iv) In the fourth figure, we set mi = 1.4 M©, m,2 — 0.16 Mq, ms = 0.80 M©, 
v = 11.76 km s _1 and take the 3 rd star to be either an MS (red) or a WD 
(green). The results can be summarized in a tabular form : 



change of variable 


result 


decrease in the value of v 


-^exchange and X m erger 

increase slightly at higher P or t 


decrease in the value of mi 


X eX change and ^merger 

increase significantly at lower P or b- 
Ionization is significant at very low 
values of 7772 and higher values of P or b 


increase in the value of 7713 


X eX change and X m erger 

increases significantly. 


nature of the incoming star 


Xexchange is slightly higher at lower P or \, when the 


(of mass 7713) 


incoming star is a WD rather than a MS. 

X merger is significantly lower when the 
incoming star is a WD rather than a MS. 



Eqn. (|23|) reveals that a oc Xa 2 n v 2 /v 2 for any fixed set of stellar parameters 
(7711,7712,7713) which becomes a oc Xa,i n /v 2 after using the expression of v c as 
in Eqn. (jlOl) . So the interaction timescale tint — 1/nav can be written as 
tint oc v/(nXai n ). As it is clear from the top-left panel of Fig. (O, that for a 
fixed set of 7711,7772,7773, the variation of X is negligible for any given <Zj n . So 
we can write U n t oc 71/77. Remember, the same relation has been used for tfi v 
in Eqns. ([2^al[2%| . 

STARLAB also gives the properties of final states, i.e. eccentricities and semi 
major axes of the final binaries for each set of inputs (mi, R\; 7772, R2', 7773, R3; a, n ). 
So for each value of P or b, in, we get a number of values of P or b. fin and e/i„. 
For each value of P or b. in, we calculated 25 percentile, median and 75 percentile 
values of P or b, fin as well as the values from analytic expressions given in Eqns. 
(TPf|) and (|15[) putting A = 0. As all these values are very close to each other, 
we use the P or b, in — Porb, fin relation corresponding to A = throughout 
(eq. Illj) . As an example, in Fig. [8l we have plotted P or b, in Vs P or b. fin 
as reported by STARLAB (scatter plots) for both merger and exchange with 
mi = 1.4 M Q , m 2 = 0.16 M Q , m 3 = 0.33 M Q , v = 11.76 km s" 1 . The 
line for 25 percentile, median, 75 percentile along with the analytic relation of 
Porb. in Porb, fin with A — (eq. Ill|) are shown in the same plot. 

In each panel of Figs. ([9]) and (flOl) . we plot P or b.in of the initial binary 
(comprising of stars mi and 7712) along the top x-axis and P or b,fin along the 
bottom x-axis. The left y-axis gives the final eccentricities while the right y- 
axis gives the time scales of the interactions. Purple points in the left panels 
are for exchange interactions while the green points in the right panels are for 
the merger interactions. Interaction time scales are plotted with black +s. The 
vertical orange lines form the boundaries of the allowed orbital period regions 
where interaction time scales are less than 10 10 years. It is clear from the scatter 
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plots (Figs, f^l ITU|) that the final binaries will most probably have e > 0.1 if they 
undergo either exchange or merger events. The observed binary pulsars with 
e > 0.1 in Ter 5 can be found in Table [5] where the companion masses (m c ) 
are for inclination angle i = 60°. These binaries are also shown in all panels 
of Figs. 191 and ITU1 (red in color, named in lower panels of Fig. [TU|) . In Fig [HI 
mi = 1.4 M , m 3 = 0.33 M Q and we vary m 2 from 0.024 M Q to 0.16 M Q and 
0.40 Mq. We have also performed the similar kind of numerical simulations by 
taking 7773 = 0.50 M Q and found that the outcomes of the simulations i.e. the 
appearance of the scatter plots (Fig [TU]) does not change much in comparison 
to the change with the variation in the value of m 2 (for fixed 7713). 

In Section |3~2TTI we have seen that PSR B1638+36B (in M13) has such a 
combination of P or b — e that it can not be explained by fly-by interactions. 
Moreover, the scatter plots (Figs. [9l [TO]) also reveal that its eccentricity can not 
be caused by exchange or mergers as these interactions imparts eccentricities 
larger than 0.01 whereas the upper limit of the eccentricity of PSR B1638+36B is 
0.001. One possibility is that the true eccentricity of this system is significantly 
lower than the upper limit so that it falls in the tji y < 10 10 yrs region. Another 
possibility is that it could have undergone eccentricity pumping of the inner orbit 
due to the presence of a third star in a wide outer orbit, i.e. in a hierarchical 
triplet like the PSR B1620-26 in M4 system [7J [7J [7J [75] . 

If an observed eccentric binary pulsar lies in the region where time scale 
for a particular interaction is greater than 10 10 years, then that interaction can 
not be responsible for its eccentricity. Moreover, an exchange interaction to be 
the origin of the eccentricity of a particular pulsar, we should have m c ~ 7773 
and for merger m c ~ m 2 + 777,3, where m c is the companion mass. Considering 
these two facts, what we can conclude about possible dynamic history about the 
eccentric pulsars in Ter 5 which are summarized below where "E" means possi- 
ble exchange and "M" means possible merger, "Ne" means exchange scenario 
satisfies the mass condition, but should be excluded on the basis of time-scale 
condition, "Nm" means merger scenario satisfies the mass condition, but should 
be excluded on the basis of time-scale condition. Values of m c s are given below 
the pulsar names. All masses are in the unit of M@. It seems that 7773 = 0.33 Af© 
is a better choice than 7773 = 0.50 Mq as the earlier one can explain the origin 
of all six eccentric pulsars in Ter 5. 
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Figure 7: Variation of X for exchange (+), merger (a) and ionization (o, whenever signifi- 
cant) with P or i, i„. (i) The first figure is for mi = 1.4 Mq, ni2 = 0.16 M©, rrt3 = 0.33 Mq 
and two values of v e.g 13.16 km s — 1 (red) and 3.27 km s _1 (green), (ii) The second figure is 
for mi = 1.4 Mq, 7713 = 0.33 Mq, v = 11.76 km s _1 and different values of m,2 e.g 0.40 Mq 
(red), 0.16 Mq (green) and 0.024 Mq (blue), (iii) The third figure is for mi = 1.4 Mq, rri2 = 
0.16 M , v = 11.76 km s" 1 and different values of m 3 e.g. 0.33 Mq (red), 0.80 Mq (green), 
(iv) The fourth figure is for mi = 1.4 M , mi = 0.16 Mq, 7713 = 0.80 Mq, v = 11.76 km s -1 
and the 3 rd star is either a MS (red) or a WD (green). 
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Figure 8: P a . 



Vs P, 



orb, fin 



plots. Right hand side (purple points) diagram is for ex- 



changes and the left hand side diagram (green points) is for mergers. The stellar parameters 
are as follows : mi = 1.4 Mq (NS), m 3 = 0.33 Mq (MS), m 2 = 0.16 Mq (MS) and 
v = 11.76 km s _1 . 
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Figure 9: Time scales (denoted by black '+') and final eccentricity distributions (scatter-plot 
of points) with initial and final orbital periods (A = in Eqn. |1H for exchange (purple points) 
and merger (green points) interactions with different stellar parameters. We plot P or i, i„ along 
the top x-axis and P or bji n along the bottom x-axis. The left y-axis gives the final eccentricities 
while the right y-axis gives the time scales of interactions. The vertical orange lines form the 
boundaries of the allowed orbital period regions where interaction time scales are less than 
10 10 yrs. The stellar parameters are as follows : mi = 1.4 Mq, = 0.33 Mq and m,2 varies 
from 0.024 Mq to 0.16 Mq and then to 0.40 3w . 
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Fi gUX6 10; Plots similar to those in Fig[9]with only difference that here 1713 = 0.50 Mq. 
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We should remember that these conditions on stellar masses are mainly 
indicative and need not to be satisfied very accurately and we allow a deviation 
of ±0.1 Mq from the exact mass conditions because: 1) a slightly different choice 
of stellar masses can give the same or similar output of the simulation and 2) the 
companion masses are not known exactly in most cases; the masses are obtained 
from the mass function in terms of m c sin i and with the assumption that the 
orbital inclination angle is 60° (Table [2]). 

3.3 Ionization 

The condition for ionization to occur is v/v c > 1 (where v c is the critical velocity 
giving three unbound stars at zero energy defined in Eq. I10[) . Fig. 1111 shows 
plots of v/v c with P or b- It is clear that ionization starts only at high values of 
Porb,in] an d a lower value of mi or a higher value of 7713 or a higher value of v 
facilitate ionization for obvious reasons of available kinetic and binding energies. 



Fi gure 11: Variation of v/v c with P or b i n for different parameters. The left panel shows the 
variation of v/v c with P rb,in for different values of m2, the middle panel shows the variation 
of v/v c with P or (, i n for different values of and the right panel shows the variation of v/v c 
with P or (, ; n for different values of v. 

STARLAB runs agree with the expectation above. As an example, keeping 
mi = 1.4 Mq and — 0.33 Mq fixed, the minimum value of P or b,in (at 
which ionization starts) increases from ~ 200 days to ~ 3000 days when we 
increase m-i from 0.024 Mq to 0.16 Mq. For vn,2 = 0.40 Mq, ionization occurs 
when P orb > 10000. 

There are three disk binaries with P orb > 1000 days such as (i) PSR B0820+20 
with P orb = 1232.40400 days, e = 0.0118689, m c = 0.23 M , characteristic 
age t c = 1.31 x 10 8 years (ii) PSR B1259-63 with P orb = 1236.72404 days, 
e = 0.8698869, m c = 4.14 M , t c = 13.32 x 10 5 years and (in) PSR J1638-4725 
with P orb = 1940.9 days, e = 0.9550000, m c = 8.078 M , t c = 2.53 x 10 6 
years. There are no GC pulsar with P or b > 1000 days. As GCs are very old 
systems, primordial binaries having one member as massive as the companions 
of PSR B1259-63 or PSR J1638-4725 would not exist still today - the massive 
companions would have evolved by now and depending upon the nature of its 
evolution, a low mass WD as a binary companion may be the end product. For 
significantly low values of the mass of the WD, the resultant binary would have 
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been ionized as have seen from FigQT]that if m 2 ^ 0.16 Mq, ionization can be 
effective for P or b ~ 1000 days. But in case of a comparatively massive WD com- 
panion, the binary will survive ionization. Although such long orbital period 
binaries can not originate through tidal capture process, there is a possibility 
of formation of such wide binaries through exchange or merger interactions (see 
Figs. l^landfTU)). Moreover, there are 12 pulsars in the galactic disk having P or b 
in the range of 100 — 1000 days, one of them has m c = 15.82 M , one has 
m c — 0.78 M Q and the other 10 have m c in the range of 0.14 — 0.38 M Q ; all 
of these have characteristic ages 10 9 — 10 10 years. But there are only 3 pulsars 
in GCs having P or b in the range of 100 — 1000 days, having m c s in the range 
of 0.21 — 0.35 Mq. As our study has shown that ionization to be operative at 
Porb < 1000 days, m 2 should be less than 0.07 M Q (which is very unlikely), 
it seems to be surprising that the number of binary pulsars having P or b in the 
range of 100 — 1000 days are significantly less than the number of these objects 
in the disk and there is no binaries with P or b > 1000 days in GCs at all. It is 
possible that such long period binaries in GCs have been misidentified as iso- 
lated pulsars as many globular clusters pulsars (e.g. those in 47 Tuc) are very 
weak, usually weaker than the detection threshold and one can see them very 
rarely only when their intensity is boosted by scintillation (M. Kramer, private 
communication). These occasions may be too rare and/or too short to identify 
them as binaries, because the larger the value of P or b, the less is the effect of the 
binary motion in the observed spin period and only observations over a longer 
time span can reveal the binary nature. 

3.4 Comments on Zero Eccentricity Pulsars 

The nearly circular binaries in GCs (group III) have P or b in the range of 
0.106 — 4.0 days. These systems lie in the e — P or b plane in such a region that 
stellar encounters are favorable. But it is clear that they have not encountered 
any kind of eccentricity imparting stellar encounters e.g. fly- by or exchange or 
merger. Dynamical formation of ultra-compact binaries involving intermediate 
mass main sequence stars in the early life of the GCs can be the origin of these 
pulsars [32] • These companions must have been massive enough (beyond the 
present day cluster turn-off mass of 0.8 Mq) so that the initial mass transfer 
became dynamically unstable, leading to common envelope evolution and sub- 
sequent orbital decay and circularization. Alternately, the present day red giant 
and NS collisions may lead to a prompt disruption of the red giant envelope 
when the system ends up as an eccentric NS-WD binary [59] . These binaries 
could decay to the group III pulsars by gravitational radiation if they had short 
Porb < 0.2 d. In any one of the scenarios mentioned above, the resultant binaries 
have to be young enough not to face any stellar encounter. In this connection, it 
will be interesting to determine the ages of these binaries by performing optical 
and/or infrared observations of their WD companions to determine WD ages 
from their temperatures. Secondly, these binary pulsars may in fact lie outside 
cluster cores with a higher viq / 714 so that the timescales for stellar interactions 
are actually > 10 10 years (even though in projection many of them appear to 
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be inside the cluster cores). 

Here we wish to mention that the positions of the binary pulsars in GCs 
with respect to the GC center in the projected image is given in Table [5J There 
are a number of binary pulsars for which position offsets are not known. Also 
the binaries which appear to be inside the cluster core in the projected image 
can actually lie outside the core. The observed values of spin period derivatives 
(P s ) of the GC pulsars can give some indication of the pulsar position and its 
environment since negative or positive P s indicate positions in the back half or 
in the front half of the cluster respectively [39]; however, the actual positions 
can not be determined with the knowledge of P s only. Dispersion measure (DM) 
can also be a tool to measure pulsar positions with respect to cluster center if 
DM variations were dominated by intra-cluster gas [76]. However, it has been 
shown that for Ter 5, DM is mainly due to the interstellar medium [76], so 
DM values alone can not help to determine pulsar positions with respect to the 
center for this cluster. 



Table 1: Globular cluster parameters taken from Webbink (1985) [77] . 



group 


GC 


d (from sun) 


via 


ri4 


v w /n 4 






(kpc) 


(10 km/sec) 


(10 4 pc" 3 ) 


(10" 3 km s- 1 pc 3 ) 


1 


Ter 5 


10.3 


1.18 


479.77 


0.0024 


2 


NGC 6440 


8.4 


1.30 


102.10 


0.013 




M 30 


8.0 


0.52 


36.39 


0.014 




NGC 1851 


12.1 


0.98 


60.40 


0.016 




M 62 


6.9 


0.52 


32.36 


0.016 




M 15 


10.3 


0.86 


36.65 


0.023 


3 


NGC 6441 


11.2 


1.44 


31.70 


0.046 




NGC 6544 


2.7 


0.59 


9.75 


0.060 




47 Tuc 


4.5 


1.32 


31.99 


0.062 


4 


M 28 


5.6 


1.06 


9.59 


0.110 




NGC 6342 


8.6 


0.45 


3.49 


0.130 




NGC 6752 


4.0 


0.78 


5.97 


0.131 




NGC 6760 


7.4 


0.58 


4.17 


0.138 




NGC 6539 


8.4 


0.45 


3.09 


0.146 




NGC 6397 


2.3 


0.48 


3.16 


0.152 


5 


M 4 


2.2 


0.51 


1.63 


0.32 




M 5 


7.5 


0.84 


2.57 


0.328 




M 3 


10.4 


0.82 


1.44 


0.568 




M 22 


3.2 


0.83 


1.96 


0.423 


6 


M 71 


4.0 


0.33 


0.35 


0.934 




M 13 


30.4 


0.78 


0.72 


1.082 




NGC 6749 


7.9 


0.39 


0.23 


1.713 




M 53 


17.8 


0.65 


0.30 


2.167 
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Table 2: Parameters for 73 binary pulsars with known orbital solutions in 23 
globular clusters. Pulsars marked with * are not used in the analysis. Pul- 
sar offsets from the cluster cores are in the units of the core radius r c (taken 
from Freire's webpage at http://www.naic.edu/~pfreire/GCpsr.html and S. 
Ransom's webpage at |http://www.cv.nrao.edu/~sransom/^ In Freire's web- 
page, eccentricities of 32 pulsars are listed as zero, which we have replaced 
by an arbitrarily small value 0.0000003 (see text for details). Last column 
shows the median values of companion masses i.e. for inclination angle i — 60 
degree (again from Freire's webpage) taking pulsar masses as 1.35 M . For 
some pulsars, better mass measurements are available from the measurement 
of rates of periastron advance or from the spectral analysis of the optical coun- 
terpart. These pulsars are listed here : (i) 47Tuc H has m p = 1.44 Mq, m c = 
0.17 Mq [78] (ii) NGC 1851A has m p = 1.35 Mq, m c = 1.103 M [79] 
(hi) M5B has m p = 2.08 M®, m c = 0.21 Mq [50] (iv) Ter 5 I has m p = 
1.87 M Q , m c = 0.30 M Q [81] (v) Ter 5 J has m p = 1.73 M , m c = 0.47 M 
[81] (vi) NGC 6440B has m p = 2.74 M , m c = 0.18 M [82] (vii) NGC 
6441A has m p = 1.26 M , m c = 0.67 M [82] (viii) PSR J1911-5958A has 
m p = 1.34 Mq, m c = 0.175 Mq (using spectral analysis of the optical counter- 
part [53]) (ix) M15 C has m p = 1.358 Mq, m c = 1.354 M [51]. 
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rc) 


(ms) 


sec/sec) 


(cm - 3 pc) 
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(Mq) 


1 


47Tuc 


J0024-7205E 


1.477 


3.536 


9.851 


24.23 


2.25684 


0.0003152 


0.18 


2 


47Tuc 


J0024-7204H 


1.75 


3.210 


-0.183 


24.36 


2.35770 


0.070560 


0.19 


3 


47Tuc 


J0024- 72041 


0.659 


3.485 


-4.587 


24.42 


0.22979 


< 0.0004 


0.015 


4 


47Tuc 


J0023-7203J (c) 


2.273 


2.100 


-0.979 


24.58 


0.12066 


< 0.00004 


0.024 


5 


47Tuc 


J0024-7204O 


0.136 


2.643 


3.035 


24.36 


0.13597 


< 0.00016 


0.025 


6 


47 Tuc 


J0024-7204P 


? 


3.643 


? 


24.30 


0.1472 


0.0000003 


0.02 


7 


47 Tuc 


J0024-7204Q 


2.227 


4.033 


3.402 


24.29 


1.18908 


0.000085 


0.21 


8 


47Tuc 


J0024-7204R (e) 


? 


3.480 


? 


24.40 


0.0662 


0.0000003 


0.030 


9 


47Tuc 


J0024-7204S 


0.432 


2.830 


-12.054 


24.35 


1.20172 


0.000394 


0.10 


10 


47Tuc 


J0024-7204T 


0.773 


7.588 


29.37 


24.39 


1.12618 


0.00040 


0.20 


11 


47Tuc 


J0024-7203U 


2.136 


4.343 


9.523 


24.33 


0.42911 


0.000149 


0.14 


12 


47Tuc 


J0024-7204V (e?) 


? 


4.810 


? 


24.10 


0.227 


0.0000003 


0.35 


13 


47Tuc 


J0024-7204W (e) 


0.182 


2.352 


? 


24.30 


0.1330 


0.0000003 


0.14 


14 


47Tuc 


J0024-7204Y 


? 


2.197 


? 


24.20 


0.52194 


0.0000003 


0.16 


15 


NGC1851 


J0514-4002A 


1.333 


4.990 


0.117 


52.15 


18.78518 


0.8879773 


1.10 


16 


M53 


B1310+18 


? 


33.163 


? 


24.00 


255.8 


0.01 


0.35 


17 


M3 


J1342+2822B 


0.254 


2.389 


1.858 


26.15 


1.41735 


0.0000003 


0.21 


18 


M3 


J1342+2822D 


0.418 


5.443 


? 


26.34 


128.752 


0.0753 


0.21 
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Table 2: (continued). 
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19 


M5 


B1516+02B 


0.545 


7.947 


-0.331 


29.45 


6.85845 


0.13784 


0.13 


20 


M5 


J1518+0204C (e) 


? 


2.484 


? 


29.30 


0.087 


0.0000003 


0.038 


21 


M5 


J1518+0204D 


? 


2.988 


? 


29.30 


1.22 


0.0000003 


0.20 


22 


M5 


J1518+0204E 


? 


3.182 


? 


29.30 


1.10 


0.0000003 


0.15 


23 


M4 


B1620-26 


0.924 


11.076 


-5.469 


62.86 


191.44281 


0.02531545 


0.33 


24 


M13 


B1639+36B 


? 


3.528 


? 


29.50 


1.25911 


< 0.001 


0.19 


25 


M13 


J1641+3627D 


? 


3.118 


? 


30.60 


0.591 


0.0000003 


0.18 


26 


M13 


J1641+3627E (e?) 


? 


2.487 


? 


30.30 


0.117 


0.0000003 


0.02 


27 


M62 


J1701-3006A 


1.778 


5.241 


-13.196 


115.03 


3.80595 


0.000004 


0.23 


28 


M62 


J1701-3006B (e) 


0.155 


3.594 


-34.978 


113.44 


0.14455 


< 0.00007 


0.14 


29 


M62 


J1701-3006C 


0.972 


3.806 


-3.189 


114.56 


0.21500 


< 0.00006 


0.08 


30 


M62 


J1701-3006D 


? 


3.418 


? 


114.31 


1.12 


0.0000003 


0.14 


31 


M62 


J1701-3006E (e) 


? 


3.234 


? 


113.78 


0.16 


0.0000003 


0.035 


32 


M62 


J1701-3006F 


? 


2.295 


? 


113.36 


0.20 


0.0000003 


0.02 


33 


NGC6342 


B1718-19 (c) 


46.000 


1004.04 


1.59 x 10 s 


71.00 


0.25827 


< 0.005 


0.13 


34 


NGC6397 


J1740-5340 (c) 


18.340 


3.650 


16.8 


71.80 


1.35406 


< 0.0001 


0.22 


35 


Ter5 


J1748-2446A (c) 


2.778 


11.5632 


-3.400 


242.10 


0.075646 


0.0000003 


0.10 


36 


Ter5 


J1748-2446E 


out(1.6) 


2.19780 


? 


236.84 


60.06 


0.02 


0.25 


37 


Ter5 


J1748-2446I 


? 


9.57019 


? 


238.73 


1.328 


0.428 


0.24 


38 


Ter5 


J1748-2446J 


? 


80.3379 


? 


234.35 


1.102 


0.350 


0.39 


39 


Ter5 


J1748-2446M 


in(0.48) 


3.56957 


? 


238.65 


0.4431 


0.0000003 


0.16 


40 


Ter5 


J1748-2446N 


in(0.39) 


8.66690 


? 


238.47 


0.3855 


0.000045 


0.56 


41 


Ter5 


J1748-24460 (c) 


in(0.45) 


1.67663 


? 


236.38 


0.2595 


0.0000003 


0.04 


42 


Ter5 


J1748-2446P (e) 


in(0.74) 


1.72862 


? 


238.79 


0.3626 


0.0000003 


0.44 


43 


Ter5 


J1748-2446Q 


out(1.45) 


2.812 


? 


234.50 


30.295 


0.722 


0.53 


44 


Ter5 


J1748-2446U 


? 


3.289 


? 


235.46 


3.57 


0.61 


0.46 


45 


Ter5 


J1748-2446V 


in(0.90) 


2.07251 


? 


239.11 


0.5036 


0.0000003 


0.14 


46 


Ter5 


J1748-2446W 


in(0.42) 


4.20518 


? 


239.14 


4.877 


0.015 


0.34 


47 


Ter5 


J1748-2446X 


? 


2.99926 


? 


240.03 


4.99850 


0.3024 


0.29 


48 


Ter5 


J1748-2446Y 


in(0.55) 


2.04816 


? 


239.11 


1.16443 


0.00002 


0.16 


49 


Ter5 


J1748-2446Z 


? 


2.46259 


? 


238.85 


3.48807 


0.7608 


0.25 


50 


Ter5 


J1748-2446ad (c) 


? 


1.39595 


? 


235.60 


1.09443 


0.0000003 


0.16 


51 


Ter5 


J1748-2446ac 


in(0.42) 


3.65859 


? 


238.75 


0.17073 


0.0000003 


0.019 


52 


NGC6440 


J1748-2021B 


0.530 


16.760 


-32.913 


220.92 


20.550 


0.570 


0.090 


53 


NGC6440 


J1748-2021D (c) 


4.230 


13.496 


58.678 


224.98 


0.286 


0.0000003 


0.14 


54 


NGC6440 


J1748-2021F 


0.690 


3.794 


31.240 


224.10 


9.83397 


0.0531 


0.35 


55 


NGC6441 


J1750-37A 


1.910 


111.609 


566.1 


233.82 


17.3 


0.71 


0.70 


56 


NGC6441 


J1750-3703B 


3.000 


6.074 


1.92 


234.39 


3.61 


0.0000003 


0.19 


57 


NGC6539 


B1802-07 


0.463 


23.1009 


47.0 


186.38 


2.61676 


0.21206 


0.35 


58 


NGC6544 


B1802-07 


? 


3.05945 


? 


134.0 


0.071092 


0.0000003 


0.010 
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Table 2: (continued). 



No 


GC 


Pulsar 


offset 


P s 




DM 


^orb 


e 










(in 




(io- 20 
















re) 


(ms) 


sec/sec) 


(cm _3 pc) 


(d) 




(M e ) 


59 


M28 


J1824-2452C 


? 


4.159 


? 


120.70 


8.078 


0.847 


0.30 


60 


M28 


J1824-2452D 


? 


79.832 


? 


119.50 


30.404 


0.776 


0.45 


61 


M28 


J1824-2452G 


? 


5.909 


? 


119.40 


0.1046 


0.0000003 


0.011 


62 


M28 


J1824-2452H (e) 


? 


4.629 


? 


121.50 


0.435 


0.0000003 


0.20 


63 


M28 


J1824-2452I (e) 


? 


3.93185 


? 


119.00 


0.45941 


0.0000003 


0.20 


64 


M28 


J1824-2452J 


? 


4.039 


? 


119.20 


0.0974 


0.0000003 


0.015 


65 


M28 


J1824-2452K 


? 


4.46105 


? 


119.80 


3.91034 


0.001524 


0.16 


66 


M28 


J1824-2452L 


? 


4.10011 


? 


119.00 


0.22571 


0.0000003 


0.022 


67 


M22 


J1836-2354A 


? 


3.35434 


? 


89.10 


0.20276 


0.0000003 


0.020 


68 


NGC6749 


J1905+0154A 


0.662 


3.193 


7 


193.69 


0.81255 


0.0000003 


0.090 


69 


NGC6752 


J1911-5958A 


37.588 


3.26619 


0.307 


33.68 


0.83711 


< 0.00001 


0.22 


70 


NGC6760 


J1911+0102A 


1.273 


3.61852 


-0.658 


202.68 


0.140996 


< 0.00013 


0.020 


71 


M71 


J1953+1846A (c) 


? 


4.888 


? 


117.00 


0.1766 


0.0000003 


0.032 


72 


M15 


B2127+11C 


13.486 


30.5293 


499.1 


67.13 


0.33528 


0.681386 


1.13 


73 


M30 


J2140-2310A (c) 


1.117 


11.0193 


-5.181 


25.06 


0.17399 


< 0.00012 


0.11 


74* 


M30 


J2140 - 2310B* 


? 


13.0 


? 


25.09 


> 0.8 


> 0.52 


> 0.02 



34 



Table 3: Binary pulsars in the galactic disk compiled from the online ATNF pul- 
sar catalog [http: / / www.atnf.csiro.au/research/pulsar/psrcat /\ Pulsars marked 
with * are not used in the analysis. Last column shows the median values of 
companion masses i.e. for inclination angle i = 60 degree taking pulsar masses 
as 1.35 M Q . 



No. 


Pulsar 


Ps 


Ps 


DM 


^orb 


c 


m c 






ms 


10~ 20 s/s 


cm -3 pc 


d 




Mq 


1 


J0034-0534 


1.87718 


0.4959 


13.76 


1.58928 


1.4E-5 


0.165 


2 


J0218+4232 


2.32309 


7.7389 


61.25 


2.02885 


2.22E-5 


0.196 


3 


J0407+1607 


25.70174 


7.9 


35.65 


669.0704 


0.0009368 


0.222 


4 


J0437-4715 


5.75745 


5.72937 


2.64 


5.74105 


1.9180E-5 


0.164 


5* 


J0610-2100* 


3.86132 


1.235 


60.67 


0.28602 


? 


0.025 


6 


J0613-0200 


3.06184 


0.9591 


38.78 


1.19851 


5.5E-6 


0.150 


7 


J0621+1002 


28.85386 


4.732 


36.60 


8.31868 


0.0024574 


0.527 


8 


B0655+64 


195.67094 


68.53 


8.77 


1.02867 


0.0000075 


0.796 


9 


J0737-3039A 


22.699378 


175.993 


48.92 


0.10225 


0.0877775 


1.559 


10 


J0737-3039B 


2773.46077 


8.92 x 10 4 


48.92 


0.10225 


0.0877775 


1.739 


11 


J0751+1807 


3.47877 


0.77860 


30.25 


0.26314 


7.1E-7 


0.150 


12 


B0820+02 


864.87280 


1.04550 x 10 4 


23.73 


1232.404 


0.0118689 


0.226 


13 


J0900-3144 


11.10965 


4.912 


75.70 


18.73764 


1.03E-5 


0.423 


14* 


J1012+5307* 


5.25575 


1.71342 


9.02 


0.60467 


? 


0.125 


15 


J1022+1001 


16.45293 


4.334 


10.25 


7.80513 


9.729E-5 


0.853 


17 


J1045-4509 


7.47422 


1.755 


58.17 


4.08353 


2.41E-5 


0.186 


18 


J1125-6014 


2.63038 


0.401 


52.95 


8.75260 


7.9E-7 


0.328 


19 


J1141-6545 


393.89783 


4.294683 x 10 5 


116.08 


0.19765 


0.171884 


1.213 


20 


J1157-5112 


43.58923 


14.3 


39.67 


3.50739 


0.000402 


1.456 


21 


J1216-6410 


3.53937 


0.162 


47.40 


4.03673 


6.8E-6 


0.182 


22 


J1232-6501 


88.28191 


81.0 


239.4 


1.86327 


0.00011 


0.166 


23 


B1257+12 


6.21853 


11.43341 


10.16 


25.262 


0.0 


0.000 


24 


B1259-63 


47.76251 


2.276484 x 10 4 


146.72 


1236.72404 


0.8698869 


4.140 


25 


J1420-5625 


34.11713 


6.8 


64.56 


40.29452 


0.003500 


0.438 


26 


J1435-6100 


9.34797 


2.45 


113.7 


1.354885 


1.05E-5 


1.079 


27 


J1439-5501 


28.63489 


14.18 


14.56 


2.11794 


4.99E-5 


1.376 


28 


J1454-5846 


45.24877 


81.7 


115.95 


12.42306 


0.001898 


1.047 


29 


J1455-3330 


7.98720 


2.428 


13.57 


76.17457 


0.0001700 


0.297 


30 


J1518+4904 


40.93499 


2.7190 


11.61 


8.63400 


0.2494845 


0.993 


31 


J1528-3146 


60.82223 


24.9 


18.16 


3.18034 


0.000213 


1.155 


32 


B1534+12 


37.90444 


242.26281 


11.61 


0.42074 


0.2736767 


1.624 


33 


J1600-3053 


3.59793 


0.94959 


52.32 


14.34846 


0.0001737 


0.240 


34 


J1603-7202 


14.84195 


1.574 


38.05 


6.30863 


9.28E-6 


0.338 


35* 


J1618-39* 


11.98731 


? 


117.5 


22.8 


? 


0.202 


36 


J1638-4725 


763.9335 


4.8 x 10 4 


552.1 


1940.9 


0.955 


8.078 


37 


J1640+2224 


3.16331 


0.28276 


18.43 


175.46066 


0.0007973 


0.290 


38 


J1643-1224 


4.62164 


1.849 


62.41 


147.01739 


0.0005058 


0.139 


39 


J1709+2313 


4.63120 


0.363 


25.35 


22.71189 


0.0000187 


0.319 


10 


J1711-4322 


102.61829 


2666.0 


191.5 


922.4708 


0.002376 


0.237 


41 


J1713+0747 


4.57014 


0.85289 


15.99 


67.82513 


0.0000749 


0.324 


42 


J1732-5049 


5.31255 


1.38 


56.84 


5.26300 


9.8E-6 


0.209 


43 


J1738+0333 


5.85009 


2.409 


33.78 


0.35479 


4.0E-6 


0.104 


44 


J1740-3052 


570.30958 


2.54969 x 10 6 


740.9 


231.02965 


0.5788720 


15.820 


45* 


J1741+1354* 


3.74715 


? 


24.0 


16.335 


? 


0.281 


46 


J1744-3922 


172.44436 


155.0 


148.1 


0.19141 


0.000 


0.097 


47 


J1745-0952 


19.37630 


9.5 


64.47 


4.94345 


1.8E-5 


0.126 


48 


J1751-2857 


3.91487 


1.126 


42.81 


110.74646 


0.0001283 


0.226 


49 


J1753-2240 


95.13781 


7.9 


158.6 


13.63757 


0.303582 


0.582 


50 


J1756-2251 


28.46159 


101.71 


121.18 


0.31963 


0.180567 


1.353 


51 


J1757-5322 


8.86996 


23(53 


30.82 


0.45331 


4.0E-6 


0.667 



Table 3: (continued). 



No. 


Pulsar 


Ps 


Ps 


DM 


^orb 


e 


m c 






ms 


10- 20 s/s 


cm -3 pc 


d 




M @ 


52 


J1802-2124 


12.64759 


7.2 


149.6 


0.69889 


3.2E-6 


0.982 


53 


B1800-27 


334.41543 


1711.3 


165.5 


406.781 


0.000507 


0.168 


54 


J1804-2717 


9.34303 


4.089 


24.67 


11.12871 


3.5E-5 


0.235 


55 


J1810-2005 


32.82224 


15.1 


240.2 


15.01202 


2.5E-5 


0.329 


56 


J1811-1736 


104.18195 


90.1 


476.0 


18.77917 


0.828011 


1.041 


57 


J1822-0848 


834.83927 


1.35 x 10 4 


186.3 


286.8303 


0.058962 


0.383 


58 


B1820-11 


279.82870 


1.37862 X 10 


428.59 


o r t *-rf i r\r\ 

357.76199 


0.794608 


0.779 


59 


J1829+2456 


41.00982 


5.25 


13.8 


1.17603 


0.1391412 


1.569 


60 


B1831-00 


520.95431 


1053.0 


88.65 


1.81110 


0.004472 


0.073 


61 


J1841+0130 


29.77277 


817.0 


125.88 


10.47163 


8.2E-5 


0.111 


62 


J1853+1303 


4.09180 


0.885 


30.57 


115.65379 


0.0000236 9 


0.281 


63 


B1855+09 


5.36210 


1.771 


13.29 


12.32717 


2.09192E-5 


0.284 


64 


J1903+0327 


2.14991 5 


1.879 


297.54 


95.17412 


0.4366784 11 


1.084 


65 


J1904+0412 


71.09490 


11.0 


185.9 


14.93426 


2.2E-4 


0.258 


66 


J1906+0746 


144.07193 


2.0280 x 10 B 


217.78 


0.16599 


0.085303 


0.976 


67 


J1909-3744 


2.947108 


1.40241 


10.39 


1.53345 


1.302E-7 


0.229 


68 


J1910+1256 


4.98358 


0.977 


38.06 


58.46674 


0.0002302 2 


0.224 


69 


J1911-1114 


3.62574 


1.416 


30.97 


2.71656 


1.9E-5 


0.140 


70 


B1913+16 


59.03000 


862.713 


168.77 


0.32300 


0.6171338 


1.056 


71 


J1918-0642 


7.64587 


2.4 


26.60 


10.91318 


2.2E-5 


0.278 


72 


J1933-6211 


3.54343 


0.37 


11.50 


12.81941 


1.2E-6 


0.382 


73 


B1953+29 


6.13317 


2.9739 


104.58 


117.34910 


0.0003303 


0.208 


74 


B1957+20 


1.60740 


1.68515 


29.12 


0.38197 


0.00000 


0.025 


75 


J2016+1948 


64.94039 


? 


34.0 


635.039 


0.00128 


0.342 


76 


J2019+2425 


3.93452 


0.70237 


17.20 


76.51163 


0.0001111 


0.364 


77 


J2033+17 


5.94895 


1.1 


25.2 


56.31 


0.00013 


0.219 


78 


J2051-0827 


4.50864 


1.2737 


20.74 


0.09911 


0.0000 


0.031 


79 


J2129-5721 


3.72635 


2.981 


31.86 


6.62549 


1.91E-5 


0.154 


80 


J2145-0750 


16.05242 


2.9757 


9.00 


6.83890 


1.928E-5 


0.503 


81 


J2229+2643 


2.97782 


0.146 


23.02 


93.01589 


0.0002556 


0.142 


82 


B2303+46 


1066.37107 


5.6909 x 10 4 


62.06 


12.33954 


0.658369 


1.431 


83 


J2317+1439 


3.44525 


0.242 


21.91 


2.45933 


0.0000005 


0.201 
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4 Eccentricity Pumping in Hierarchical Triple 
Systems 

In the previous section, we have discussed how stellar interactions in globular 
clusters lead to binary millisecond pulsars in eccentric orbits. In the absence of 
such stellar interactions i.e. in the galactic disk, millisecond pulsars are expected 
to be in circular orbits as during the accretion phase (recycling process), tidal 
coupling leads to orbit circularization. On the other hand, if the millisecond 
pulsar is the part of a hierarchical triple, consisting of inner and outer binaries, 
then as a result of Kozai oscillation, the eccentricity of the inner binary can be 
large. 

To understand this mechanism, let us assume that a star of mass too forms 
the inner binary with its companion of mass m\. Another star of mass 777,2 forms 
an outer binary with the center of mass of the inner binary. One can consider the 
inner and outer binaries to be isolated binaries of masses too, m\ and Too + mi, 
to 2 respectively. Let us denote the eccentricities, semi-major axes, inclinations 
with respect to the sky plane and the arguments of the periastrons (with respect 
to their lines of nodes) of the inner and outer binaries by ei, e2, 01, 02, ii, ii, Wi 
and u>2, respectively and the mutual inclination angle between the two orbits 
4 = i\ + «2- Then the magnitudes of the angular momenta Qi and G2 of the 
inner and outer binaries will be as follow : 

g 1 =rHiGM?a 1 (l-4)\ , (24a) 

where G is the gravitational constant, Mi = (m + mi) is the total mass of the 
inner binary and rji is the symmetric mass ratio of the inner binary, given by 
r/i = too toi /Mf. The temporal evolution of the inner binary can be probed with 
the help of secular perturbation theory, applicable to Newtonian hierarchical 
triples containing point masses, while including the dominant quadrupolar order 
interactions between the two orbits. In other words, the dynamical equations 
invoked are accurate to order (ai/a 2 ) 2 , where (ai/a 2 ) is the small parameter in 
the perturbative expansion. The relevant equations providing secular temporal 
evolution for the eccentricity and the argument of periastron of the inner binary 
including the dominant order general relativistic effect that causes the periastron 
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of an isolated compact binary to advance are [85] : 
u x = ^^Ue 2 + (5cos2wi - 1)(1 



Gi 
6C 2 9 i 

— t—< 2 + e 2 (3 - 5cos2wi) 
G: 



2 



+ 



c 2 ^(l-el) 



G Mi 



-i 3/2 



Ol 



(25a) 



on (~i 

ei = ex (1 - ef ) (1 - 2 ) sin 2wi , (25b) 



ft 

where = cosi m . The quantity C 2 is 



2 



16 a 2 (l-e^y/ 2 \a 2 J 

Note that the Newtonian contributions to Eqns. (I2~5j) originate from certain 
'doubly averaged' Hamiltonian, which is derivable from the usual Hamiltonian 
for a hierarchical triple at the quadrupolar interaction order [86]. The 'doubly 
averaged' Hamiltonian, suitable for describing secular (long-term) temporal evo- 
lution of a hierarchical triple, is independent of the mean anomalies of the inner 
and outer orbits. This implies that their respective conjugate momenta and 
hence the semi-major axes, a\ and a 2 , are constants of motion. Here changes 
in a\ and e\ due to gravitational wave radiation have been neglected. First two 
terms in Eqn. (|25aj) give the rate of change of periastron due to the hierar- 
chical triple configuration (uii^t)- The third term gives the rate of change of 
periastron due to the general relativistic effect (wi, ffr ) up to 1PN order only and 
neglecting the spin-orbit coupling (see Eqn. [1]). 

It turns out that if the mutual inclination angle is fairly high and in a cer- 
tain window, the orbital eccentricities experience periodic oscillations over time- 
scales that are extremely large compared to the respective orbital periods. The 
above effect, usually referred to as the Kozai resonance, arises due to the tidal 
torquing between the two orbits [57]- The Kozai resonance can force initially 
tiny eccentricity of the inner binary to oscillate through a maximum value, given 

1 /2 

by ei max ~ (l — | cos 2 i m ,o) where i m> o is the initial value for the mutual 
inclination angle. Due to the obvious restriction, namely |cosi m! o < (3/5) x ^ 2 1 , 
i rn fl is required to lie in the range 39° — 141°. An approximate expression for 
the period of the Kozai oscillation reads [55] : 

r K ozai=Po, 6 .i^ (-V (1-e 2 .) 3 / 2 , (27) 

where P or b,i is the orbital period of the inner binary. 

The general relativistic periastron advance of the inner binary, in princi- 
ple, can interfere with the Kozai resonance and even terminate the eccentricity 
oscillations. This is because the extra contribution d)i tgr to uix can indirectly 
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affect the evolution of e±. The following useful criterion can be used to infer the 
possibility for this not to happen |85j : 



< 



3 mi a i ( 1 



3/2- 



1/3 



(28) 



where M< is Mj/M , 

o-i = a i/Lo and Lq = 1.476625 km. 
Another constraint for 0,2/0,1 can be obtained by invoking an empirical re- 
lation needed for the Newtonian coplanar prograde orbits in hierarchical triple 
configurations [89] : 

\aij l-e-2 

We treat the above inequality to be rather conservative [85] as the inclined or- 
bits, relevant for our investigation, are expected to more stable than the coplanar 
triples of Eq. (j29|) . 



m 2 
M 



(1 + ea) 
(l-e 2 )Va 



(29) 



4.1 PSR J1903 + 0327 : Open Question 

Recently a millisecond pulsar (P s ~ 2.15 ms) PSR J1903 + 0327 in a highly 
eccentric orbit around a solar mass companion in the galactic plane has been 
discovered [43]. The discoverers almost discarded the "born- fast" scenario for 
this pulsar and strongly favored a binary origin for it. So they had to explore 
unconventional formation scenarios. They favored the idea that this pulsar is a 
member of the inner binary of a triple system and Kozai resonance is responsible 
for its high eccentricity. They have also speculated a globular cluster origin for 
the system, involving the ejection of the system from a globular cluster or the 
disruption of the cluster itself. 

Accurate radio timing analysis of the pulsar (using 'DDGR' timing model 
[90] ) reveals that the binary has the following properties - the orbital period is 
95.1741176(2) days, eccentricity is 0.436678411(12), pulsar mass is 1.74(4) M , 
companion mass is 1.05(15) M & , the longitude of periastron is 141.65779°, the 
advance of periastron is 2.46(2) x 10 -4 deg yr _1 [43]. Here the numbers in the 
parenthesis are twice the 1-tr uncertainties in the least significant digit(s). The 
companion mass was determined using the advance of periastron and the shapiro 
shape ("s") parameter, assuming the advance of periastron to be only due to the 
general relativistic effect i.e. the third term in Eqn. (j25a[) . Infrared observations 
with Gemini North telescope yielded a possible main-sequence companion star 
for the pulsar. 

The hierarchical triplet hypothesis for this system is that the pulsar and 
the companion (either a neutron star or white dwarf as the pulsar is recycled 
- implying past accretion from the red-giant phase of its companion) forms the 
inner binary and the main-sequence star is in a much wider orbit forming the 
outer binary. This means that the white-dwarf in the inner binary contributes 
in the timing analysis and the main-sequence star in the outer binary has been 
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observed in infrared observation. So, in Eqns. (|25a[) . (I25b[) . (|2"6"j) . (j2"T)) , 

|28|) and ([29]) one can use m = 1.74 M Q , ra x = 1.05 M , e x = 0.436678, 
wi = 141.658°, wi = 2.46(2) x 10~ 4 deg yr~ 4 . a\ can be calculated from the 
orbital period (~ 95.174 days) using Kepler's law. But the problem here is 
the companion mass was found to be 1.05 Mq with an assumption that loi is 
only due to general relativistic effect (third term in Eqn. I25a[ LOi :gr ) on the 
other hand, for a hierarchical triple system all the terms in Eqn. (|25al) would 
contribute. So if this system is the part of a hierarchical triple, the value of 
the companion mass of the inner binary is likely to be different from the above 
mentioned value. 

So the timing-analysis of this system was continued. The latest value of 
the companion mass measured using only the Shapiro shape ("s") and range 
parameter ("r") is 1.03 ± 0.04 Mq. This matches with the earlier estimate 
using s and u>i supporting the assumption that the ui\ is purely relativistic [4"5] . 
However, because of the lower precision in the measurement of r, it is impossible 
to exclude small contributions to uj\ at present. Using Eqn. I25al we found that 
for preferable combination of a = 02/01, 7712, &i the contribution form the first 
two terms in Eqn. I25al can be negligible in comparison to the third term. As an 
example, in Fig. [T2]we plot the variation of uj\ y ht and LOi, gr with 1712 and for a 
fixed a — 600 (for such combinations of , &i and a that both the inequalities 
(|2"5|) and (|2"9"|) are satisfied); red points are for u)\,ht and green points are for 
^>\,gr- ^i,gr is independent of mi and ei (see the third term of Eqn. I25al) and 
matches with the observed value [43]. We find that u)\.ht can be significantly 
less than wi iSr for sufficiently small values of e2. In that case the companion 
mass measured using s and r can match with the measurement done with s and 
cL>i assuming wi = cL>i tflr . But the precision of r is still improving significantly 
with continued timing and it might be possible in future to confirm whether 
there is even any small deviation of 6j\ from wi iffr . 

In principle, t\ measurement can also help to test the validity of hierar- 
chical triple scenario. For a hierarchical triple system, the observed k\ can be 
the combined effect of three terms like e.\ t0 bs — &i,ht + &i,gw + ei.geoj where 
ei m is the rate of change of eccentricity because of Kozai oscillation in a hi- 
erarchical triple system given by Eqn. (|25bp . k\^ gw is the rate of change of 
eccentricity because of gravitational wave radiation given by Eqn. (|17[) and 
ei.geo is due to geodetic precession of the pulsar spin axis (see Eqns. 8.58, 
8.82 and 8.91 of 0). It is clear from Eqn. (fl"7| that ei iffu , is always nega- 
tive and wc found it to be -1.45975808 x 10" 24 s _1 for the inner binary as 
mentioned above. Eqn. (|25b[) shows that ei^t c& n be positive or negative de- 
pending on whether sin2a>i is positive or negative and for uj\ = 141.65779°, 
sin2w 1 = —0.97312 making ei^t negative. So the magnitude of e-i^t should 
be such that (-1.45975808 x 10~ 24 s" 1 + e 1M + ei. geo ) = e Mbs . The problem 
here is that although we can calculate ei^t using Eqn. (|25b[) [4"4] using the 
parameters of the inner binary as mentioned earlier and varying i m , 02 , &i and 
7712 over realistic ranges satisfying conditions (f28|) and (|29|) . we can not calculate 
ei,geo as the orientation of the spin axis of the pulsar with respect to the total 
orbital angular momentum is not known, which is needed to calculate e\. geo 
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Fi gure 12: Variation of il)^ and 6j gr with m<i and e2 for a = 600. Red points are for oj^ 
and green points are for independent of m2 and e2 (see the third term of Eqn. 

I25ali and matches with the observed value 43 . Points are plotted only for such combinations 
of rri2 and e2 that inequalities 1281 1 and I I29II are satisfied. 

(see Eqns. 8.58, 8.82 and 8.91 of [6). Moreover, the presently determined value 
of ei from timing analysis is (8 ± 15) x 1CP 17 s _1 which is merely the mea- 
surement uncertainty (unpublished, provided by Paulo Freire through private 
communication) . So this approach can not help us to test the hierarchical triple 
hypothesis at present. 

On the other hand, the observed value of the spin period derivative P s can 
be used to test the hierarchical triple scenario. A rough estimate for P s can be 
obtained by treating the inner binary as a single object of mass mo + mi which 
forms a wide binary with stellar mass companion of mass m-i [91] . Further, 
we let the projected semi-major axis of this new binary to be 600 times that 
of the inner binary and choose a value for e 2 consistent with inequalities (|28p 
and (|29|) . The associated line of sight acceleration leads to an estimate of P s 
of the order of 10~ 17 ss -1 [6] and this is much higher than the reported value 
P s ~ 10 _20 ss _1 [43 . This is fact seems to be a clear indication that this system 
is not the part of a hierarchical triple. 

There is an alternative hierarchical triple scenario [131 UJ5], such that the 
original inner companion (WD) of the pulsar has been evaporated by now and 
we see only the outer binary at present with mo — 1.74 M Q , 7712 — 1.051 M@, 
io 2 = 141.65°, e 2 ~ 0.436678, uj 2 = 2.46(2) x 10" 4 deg yr" 1 and P orbi 2 = 95.174 



d(o/dt HT 
do)/dt GR 
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days. This scenario is difficult to test as we don't have any information about 
the parameters of the extinguished inner binary like mi, e\, Ux and oi. In this 
case, the agreement between the measurement of the present day companion 
mass (777.2) using s and 0J2 with the same using s and r is expected as u>2 is only 
due to general relativistic effect (third term in Eqn. I25a[) because there is only 
one binary at present. 




0.01 0.1 1 10 100 1000 



Porb(d) 

Fi gure 1 3 1 Fly-by isochrones for tfly = 10 years in the eccentricity - orbital period plane 
for different values of the GC parameter vio/m. The region to the lower right of each of these 
lines have shorter timescales. This indicates that PSR J1903+0327 could have originated only 
in a rather special GC environment with characteristics belonging to the upper-left portion 
of the isochrone with uio/714 = 0.15 (the unit of uio/ri4 is 10 — 3 km s — 1 pc 3 as mentioned 
earlier). 

Let us now comment on the possible globular cluster origin for the system, 
motivated by the fact that in globular clusters, stellar interactions can impart 
eccentricity to millisecond pulsar binaries and indeed ~ 20% of the millisecond 
binaries in globular clusters have eccentricity > 0.1. Using the formalism dis- 
cussed in the earlier section for fly-by interaction fsee !22bp . we find that the PSR 
J1903+0327 like systems could have been produced in globular cluster systems 
up to a mean 1)10/774 = 0.15 but not in any globular cluster with a higher value 
of this parameter. Here we consider maximum allowed timescale for stellar in- 
teraction to be 10 9 years as the characteristic age of the pulsar is 1.8 x 10 9 years 
which is usually an upper limit of the true age of a radio pulsar [43 . Thus the 
putative origin of PSR J1903+0327 inside a globular cluster could have taken 
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place only in a restricted set of globular clusters. This further restricts the envi- 
ronment for production of such pulsars and therefore the overall probability for 
origin of this pulsar in a globular cluster must be smaller than that estimated 
by the dicoveres [13] ■ Here it is noteworthy to mention that the other types of 
eccentricity inducing interactions e.g. exchanges and mergers also scale with the 
parameter vio/n^ - so these interactions also can lead to PSR J1903+0327 like 
systems only for a restricted set of globular clusters. So the eccentricity of this 
system still remains to be a puzzle to us. 

5 Summary 

In this article we discuss how orbital parameters of binary radio pulsars reveal 
the history of their formation and evolution including dynamic interactions with 
other objects, the physics of the matter at extreme densities as inside the pul- 
sar as well as the properties of the pulsar magnetosphcrc. Firstly we try to 
understand the importance of studying DNSs. If at least one of the NS of a 
DNS is a radio pulsar, we can determine the moment of inertia of the pulsar 
through accurate radio timing analysis by measuring the extra advancement 
of the periastron of the orbit due to the spin-orbit coupling over the standard 
post-Newtonian one. Then by matching this observed value of the moment of 
inertia with the theoretically predicted value, one can constrain the dense mat- 
ter EoS. One necessary condition for this to be done is that the pulsar must be 
rapidly spinning and its companion NS must be much slower; the double pulsar 
system PSR J0737-3039(A/B) seems to be the most potential candidate for this 
purpose till date as we know that PSR J0737-3039A has P s = 22.7 ms whereas 
PSR J0737-3039B has P s = 2.77 s. Moreover the masses of both the NSs are 
well determined. We also checked how good the system really is for this purpose 
by studying the effect of various parameters on the spin-orbit coupling term as a 
larger value of this term will facilitate the moment of inertia measurement. We 
found that this term has larger effect if the following conditions are satisfied - (i) 
short spin period (ii) short orbital period, (iii) high eccentricity. Unfortunately, 
PSR J0737-3039 lacks conditions (iii) and even (i) to some extent. So it will be 
nice if we can find better systems satisfying all the above conditions. Globular 
clusters are potential sources for favorable DNS systems in the above sense as 
DNSs can be formed by exchange collisions due to the high stellar density and 
the abundance of millisecond pulsars in globular clusters. On the contrary, only 
one DNS has been found in a GC so far which might be a result of selection 
effects as GCs are usually at larger distances in comparison to the disk pulsars. 
The resulting binaries are more probable to be highly eccentric with a millisec- 
ond pulsar facilitating moment of inertia estimation. We have also discussed 
how the study of the eclipse of the radio beam of PSR J0737-3039A by the 
magnetosphere of PSR J0737-3039B can be used as a tool to understand the 
magnetospheric properties of a pulsar in a better extent. Some suggestive works 
about the formation of PSR J0737-3039A/B have been mentioned too. 

We explain the presently observed orbital eccentricity and period data of GC 
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binary pulsars by stellar interaction scenarios of fly-bys, exchanges and mergers 
with held stars. Binaries with e > 0.1 are most probably the result of exchange 
or merger events whereas binaries with 0.01 > e > 0.00001 are products of fly-by 
of single stars. A number of wide orbit intermediate eccentricity pulsars seen in 
the galactic disk are absent in the GC sample because they have been kicked up 
to relatively high eccentricities by passing stars in the dense stellar environments 
in GCs. In some GCs such as Ter 5, the stellar densities are so high, and the 
velocity dispersion so modest that the interaction timescale for exchange and 
fly-by interactions is relatively short. In such GCs a typical binary system 
may undergo multiple interactions. If the original binary contains a spun-up 
millisecond pulsar in a relatively "soft" binary, then the exchange interaction 
may even produce a single millisecond pulsar in the cluster. This may explain 
the higher incidence of isolated millisecond pulsars in GCs compared to that in 
the galactic disk [33] ; in GCs 52 out of 59 isolated pulsars are millisecond pulsars 
whereas in the disk 23 out of 1565 isolated pulsars are millisecond pulsars. In 
addition, exchange interactions, as we have seen, can lead to highly eccentric 
orbits and the system can be ejected from the cluster core. If the last encounter 
took place not too long ago, the system can be at a relatively large offset from the 
cluster core, albeit being still spatially co-located with the GC. We have found 
that the eccentricity of PSR B1638+36B (in M13) is difficult to explain by any 
type of interactions - fly-by, exchange or merger. So we propose that either its 
true eccentricity is lower than the presently known upper limit or it is a member 
of a hierarchical triple. We have also considered the effects of collision induced 
ionization on the present day distribution of orbital parameters of radio pulsars 
in GCs. In the galactic disk we find 12 pulsars with 100 < P or b < 1000 days 
(although few of them have more massive companions than the pulsar binaries 
in the GCs), while there are only 3 pulsars in this range of orbital period. 
Although it is tempting to speculate that these systems are missing from the 
present day GCs because they have been ionized in the past, we find that the 
ionization probability becomes substantial in this orbital period range only for 
very small values of companion masses. Observational selection effects can be 
the cause of this discrepancy between the number of large orbital period pulsars 
in the disk in comparison to that in GCs. Many galactic disk binary pulsars are 
seen in the region of e — P or b plane predicted by Phinney due to the fluctuation 
dissipation of convective eddies and the resultant orbital eccentricities that are 
induced. These pulsars are missing from the GC sample. There is no reason 
not to expect these systems to form in the GCs (although due to the lower 
metallicity of the stellar companions in GCs, they are expected to be in slightly 
different region in the e — P or b plane) . This can be explained by the substantial 
probability of them being knocked out of their original phase space due to stellar 
interaction in GCs. Indeed there are some wide binaries in the present day GCs 
with moderately high eccentricities (e.g. 0.01 < e < 0.1 and 60 < P or b < 256 d) 
which could have arisen out of fly-bys or exchanges from progenitor binaries 
with "relic" eccentricities e ~ 10~ 4 (merger in unlikely because of moderately 
low values of m c s). The circularity of group III pulsars in GCs are thought to 
be due to either their special formation channels and young age or due to their 
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locations at the outskirts of the GCs. 

PSR J1903+0327 is the only one eccentric millisecond pulsar in the galactic 
disk discovered so far. Its millisecond spin-period suggests the spin-up scenario 
through accretion in the past. As such recycled pulsars are expected to be 
in circular orbits due to tidal-coupling induced circularization, this system has 
created lots of interest among pulsar researchers. As it is well known that 
in globular clusters, stellar interactions can impart eccentricity to an initially 
circular orbit, the possibility of past GC association of this system have been 
studied and almost rejected. Moreover, the possibility of this system being a 
member of a hierarchical triple system also seems to be less unlikely. So this 
system remains an open question for us. 

We hope in the future with more advanced technologies like SKA, many 
more new interesting types of pulsars will be discovered and that the exciting 
and challenging situation will ultimately lead to a better understanding of the 
physics of pulsars. 
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